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Abstract 



The shapes of the inclusive photon spectra in the processes J/ip — * 7X and 
T — * 7A have been analysed using all available experimental data. Relativistic, 
higher order QCD and gluon mass corrections were taken into account in the fitted 
functions. Only on including the gluon mass corrections, were consistent and ac- 
ceptable fits obtained. Values of 0.721±g;g^| GeV and 1.18±g;§| GeV were found for 
the effective gluon masses (corresponding to Born level diagrams) for the 3/ 1(1 and T 
respectively. The width ratios T(V — > hadrons)/T(V — > 7 + hadrons) V = T 
were used to determine a 5 (1.5GeV) and a s (4.9GeV). Values consistent with the 
current world average a s were obtained only when gluon mass correction factors, 
calculated using the fitted values of the effective gluon mass, were applied. A gluon 
mass ~ 1 GeV, as suggested with these results, is consistent with previous analyt- 
ical theoretical calculations and independent phenomenological estimates, as well 
as with a recent, more accurate, lattice calculation of the gluon propagator in the 
infra-red region. 
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1 Introduction 



As suggested by the inventors of QCD [[[], the colour symmetry of the theory is, 
conventionally, assumed to be unbroken, so that, theoretically @, the gluon is supposed 
to have a vanishing mass. It was also conjectured, by the same authors, that the resulting 
infra-red divergences of the theory at large distances ('infra-red slavery') might explain the 
confinement of quarks. As is also well known, in the contrary case that gluons are massive, 
there is a possible breakdown of renormalisability as well as violation of unitarity at high 
energy by certain tree level amplitudes. These problems are common to all non-abelian 
gauge theories with massive vector mesons |4], |5| . 

These problems may be solved, as in the Standard Electroweak Model, by the intro- 
duction, also for the strong interaction, of spontaneous symmetry breaking and the Higgs 
mechanism ||. Since, however, there is no experimental evidence for the existence of a 
Higgs boson for the strong interaction, or for electrically charged gluons, which are also 
predicted by some of these 'broken colour' theories, it is still generally supposed, in spite 
of the infra-red divergent nature of such a theory, that the QCD colour symmetry remains 
unbroken. 

A possible way out of this dilemma (infra-red divergences if the gluon mass is zero, 
breakdown of renormalisabilty and unitarity if it is not) was suggested by Cornwall [^, ^ 
who suggested that non-vanishing gluon mass might be dynamically generated in a theory 
in which the colour gauge symmetry remained unbroken. Other authors || |TU| pointed 
out that a gauge invariant, renormalisable and infra-red finite, version of QCD with 
massive gluons is possible, provided that a suitable four- vertex Faddeev-Popov ghost field 
is introduced into the theory. 

The aim of the present paper is not to pursue further these theoretical considerations^, 
but rather to seek direct experimental evidence on the mass of the gluon. The processes 
considered, the radiative decays of ground state vector heavy quarkonia into a photon and 
light hadrons, are particularly well adapted to such a study, as the observed final state 
results from the hadronisation of a pure two gluon final state at the lowest order in per- 
turbative QCD (pQCD). These are the 'golden' physical processes for the determination 
of the gluon mass, that may be compared to the neutral kaon system for the study of 
CP violation or Tritium /3-decay for the direct determination of the mass of the electron 
antineutrino. Indeed, the analogy between the process J/tp — > 7X and Tritium /5-decay 
is a very close one. In both cases it is the study of the end-point region of a spectrum ( 
that of the electron for Tritium /5-decay, of the photon for the radiative J/i/j decay) that 
give the mass limits on the V e or the gluon mass. The J/ip, being the lightest quarkonia is 
the most sensitive to the gluon mass, just as the Tritium /3-decay, with a very low energy 
release, gives the best direct limit on the V e mass. Indeed, as will be shown below, the 
suppression of the spectrum end point due to gluon mass effects, is much more severe in 
the case of J/ijj radiative decays than for the heavier T state. 

Already in 1980, Parisi and Petronzio [0(PP) had suggested a mass of ~ 800 MeV for 

1 The interested reader is referred to Reference for recent developments, and citations of the 
related literature 
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the gluon, on the basis of the strong suppression of the end point of the photon spectrum 
in radiative J ftp decays, as measured by the MARK II collaboration fl3| . In order to 
relate, however, in a precise way, the shape of the photon spectrum to the gluon mass, 
two other important physical effects, which also soften the shape of the photon spectrum, 
must also be properly accounted for. These are: (i) relativistic corrections and (ii) higher 
order QCD corrections. Because of the only recently available complete next-to-leading- 
order (NLO) pQCD calculation of the photon spectrum in the decays T — > 7X [111] and a 



much improved understanding of the phenomenology of relativistic corrections based on 
several recent and independent potential model calculations, the analysis presented below 
is the first to take fully into account the important effects (i) and (ii) and so confirm the 
conclusion of PP that the gluon mass is ~ 1 GeV. At the time of writing, no calculation 
yet exists in which the effects (i) and (ii), as well as that of the gluon mass are taken 
into account in a unified way, so the present analysis is inevitably a phenomenological one 
where the three different types of corrections are assumed, loosely speaking, to 'factorise'. 
Since, however, it is clear that the gluon mass effects are, by far, the most important, no 
large uncertainity on the results obtained are expected to result from this approximation. 

The results presented below also confirm the conclusions of two previous, closely re- 
lated, papers written by M.Consoli and the present author |[J, [RJ(MCJHF1,MCJHF2). 



Some brief comments are made here on these papers: some more detailed remarks are 
made in Section 8 below. 

In MCJHF1 effective gluon masses, m g , determined from fits to J/ip — > 7X and 
T — > 7 A were used, in conjunction with gluon mass correction factors calculated by 
PP (or re-calculated using pure phase-space considerations) to derive a large number of 
a s values from different charmonium and bottomonium branching ratios. Agreement with 
the expected pQCD evolution of a s from the scale Mj/^/2 to My/2 was only obtained 
when the gluon mass corrections were applied. Also, only in this case, was good agreement 
found between the derived values of a s and those obtained from deep inlastic scattering 



experiments. In this paper only the Photiadis |17| higher order (HO) QCD correction 
(that is only applicable in the end-point region of the photon spectrum, and does not 
include real gluon radiation effects) was used, and relativistic corrections were completely 
neglected. 

The second paper, MCJHF2, made essential use of the recently proposed non-relativistic 
quantum chromodynamics (NRQCD) formalism of Bodwin, Braaten and Lepage (BBL) ||18|| , 
in which both non-relativistic and HO QCD corrections (but not gluon mass effects) were 
treated in a rigorous way, order-by-order in perturbation theory, using an Operator Prod- 
uct Expansion. As suggested by BBL, the values of a s , and the leading relativistic cor- 
rection parameter r ~ (v 2 ) were treated as free parameters in fits to various charmonium 
and bottomonium decay widths. Similar fits were also performed to the inclusive photon 
spectra in J/ip and T decays. No consistent values of r and a s were found in the absence 
of gluon mass corrections. When the latter were included, consistent values of a s similar 
to those found in MCJHF1 were obtained. However, in this case, the values of r were 
found to be much smaller in absolute value than the expectations from potential model 
calculations, and even (as discussed further in Section 8 below) of the wrong sign. The 
conclusion concerning the inability of the NRQCD formalism to describe the experimental 
data, in the absence of gluon mass corrections, was not however affected by the incorrect 
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Reference 




T 


Bradley [|9[ 


0.44 


0.069 


Eichten el al. 


po| 


0.20 


0.096 


Bergstrom el al. I |31[ 
V(r) = 0.2r - ^ 
V(r) = \uu) 2 r 2 

" V 1 2 1 


0.47 
0.21 
0.34 


— 


Bergstrom el al. II |32| 
V(r) = 0.2r - ^ 
V(r) = 0.163r- ^ 


0.47 

— 


— 
0.47 


Beyer e/ a/. [28] 


0.21 


0.18 


Chiang el al. p5| 


0.21 


0.078 


Chao el al. p6| 


0.26 


0.13 


Schuler 27 
l/(r) = Ar 1 - 
i/ = -0.1 
v = 0.0 
v = 0.3 




0.36 
0.32 
0.25 


0.075 
0.066 
0.048 



Table 1: Estimations of < v 2 > for the J/ip and the T 



treatment of relativistic corrections. 

The structure of this paper is as follows. Sections 2, 3 and 4 are devoted to descriptions 
of the implementation of relativistic, HO QCD and gluon mass corrections respectively. 
Fits to the experimental data on J/ip — > 7X and T — > 7X to obtain, in each case, 
the corresponding effective gluon mass m g , are described in Sections 5 and 6. Section 
7 describes the determination of a s (m c ) and a s (mt,) from the experimental branching 
ratios T(V — > hadrons)/r(V — > 7 + hadrons), V = J/ip, T. The values of a s obtained in 
this way are unaffected by relativistic corrections. In Section 8, the effective gluon mass 
values obtained in this paper are compared with other estimates of the gluon mass in the 
literature. Finally Section 9 contains a brief summary and outlook. Details of the method 
used to simulate the effects of experimental resolution on the inclusive photon spectrum 
are given in an Appendix. 



2 Relativistic Corrections 



Relativistic corrections to the van Royen-Weisskopf formula []T9[ for the decay rate of 
a vector meson V into a charged lepton pair: 



r( v.«>^— S| #) p (2.1; 



16vra(Mv) 2 4 

were calculated by Bergstrom et al. |20|| . A relativistic correction factor, f^c, to the 
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ldT 2 

rdz 



Figure 1: Theoretical predictions for relativistic corrections to the inclusive photon spec- 
trum in a) J/ip, and b) T , radiative decays. In a), the LO QCD prediction is shown as 
the dot-dashed line. The solid line gives the prediction of Eqn(2.8) for < v 2 >= 0.28. The 
dashed line shows the prediction of Eqn(2.8) neglecting the order < v 2 > 2 term. In this 
case only the positive part of the prediction is shown. In b), the LO QCD prediction fiSSjJ 
is again shown as the dot-dashed line. The solid line gives the prediction of Eqn(2.8) for 
< v 2 >= 0.09. In both a) and b) the peak of the function in Eqn(2.8) at small values of 
z has been suppressed 
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leptonic decay width r(V — > £ + £ ) was found with the general form 

1 



f R c(V^£ + £-) = (l-^r) 2 (2.2) 



where 

d 3 p (E(p) - m Q ) ip(p) 
(2vr) 3 E(p) V(0) 



(2.3) 



and E, p and mQ are the energy, momentum and mass of the bound heavy quark. if)(p) 
is the wavefunction in relative momentum space, related to the spatial wavefunction at 
the origin, ^(0), by the expression: 

'^ (0) = / Wf^ pl {2A) 

In the approximation where the valence quarks of the vector meson are considered to be 
symmetrically bound in the meson rest frame, so that E(p) = My/2, it follows that 



E{p) "My ~ ^ 

pl +0(p i ) 



v 



,2 



/r +0(p 4 ) 



2m 



2 



V 



2 



+ 0(v 4 ). (2.5) 



Here e is the 'binding energy' My — 2rriQ and v is the velocity of the heavy quarkQ. Using 
Eqn(2.5), (2.2) may be written as: 

f*c(y - z + n = {i- l -<v 2 >f + o{v & ) (2.6) 

o 

where < v 2 > is the mean value of the squared velocity, that depends on the bound state 
potential. Similar relativistic corrections were calculated for several decay processes of 
heavy quarkonia by Keung and Muzinich (KM) [pl|. The calculation was based on a 



non-relativistic reduction of the Bethe-Salpeter equation |22| for the relativistic quark- 
antiquark bound state problem. The results of KM were presented as 0(v 2 ) corrections to 
the decay rate rather than to the decay amplitude, as in Eqn(2.6) above. In the present 
paper all relativistic corrections are applied at the amplitude level so that additional 
0(v 4 ) terms are added to the results quoted by KM to 'complete the square' and obtain 
a positive definite decay rate. This correction is important only for charmonium decays 
where, because of the relatively large value of < v 2 >, the corrected decay rate becomes 
negative for both small and large values of z = 2E 7 /My, if only the 0(v 2 ) correction 
terms are retained. KM confirm the relativistic correction factor for V — > £ + £~ given in 
Eqn(2.6) and find also: 

fnc(y -> ggg) = /*c(V - 7 gg) = (1 - 2.16 < v 2 >) 2 + 0(A (2.7) 



2 In units with c = 1 
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Of particular importance for the present study, KM also give, in their Eqn(3.5), the 
relativistic correction to the inclusive photon spectrum in V — > 7gg. 'Completing the 
square' to obtain a positive definite differential decay rate yields the spectrum: 



1 dT 

r dz 



1 

Cn 
1 

Cn 



g(z) <v z > 



f {z)+g(z) <v 2 >^ 



{g{z) < v 2 >f 



(2- 



where [F2B 



C 



N 



(vr 2 - 9) 



1+ < v 2 > 



1 (9tt 2 - 68) 
4 (tt 2 -9) 



~ (tt 2 -9)(1 -4.32 < v 2 >), (2.9) 



and 



9W 



5fo(z) h{z) 



(2.10) 



3 12 

The functions fo(z) and f\{z) are reported in Eqns(15) and (16) of Reference (T^j. In 
the approximation used here, relativistic corrections are completely specified by the sin- 
gle parameter < v 2 >. Although one may hope, in the future, to determine this non- 
perturbative parameter by lattice QCD methods [23|, the only existing estimates are 
derived from potential models of the quarkonium bound state. Some of the estimates of 
< v 2 > for the J/ifj and the T, that have been given in the literature, are presented in 
Table 1. Usually in these papers the relativistic correction factor for the decay V — > £ + £~ 
is quoted. For the entries in Table 1, this is converted into a value of < v 2 > using 
Eqn(2.6). In the case of Chiang et al. |25| Eqn(2.7) is used, and for Chao et al 
ratio fnc(y — > ggg)//i?c(V — > £ + £~). Schuler [^7] calculated directly values of < v 



26 



the 
> for 

a series of different charmonium and bottomonium states as a function of the prameter v 
in a power-like potential of the form V(r) = \r u . The range of different values of < v 2 > 
presented in Table 1 is very wide: 0.20-0.47 for the J /if) and 0.048-0.47 for the T. Apart 
from the estimates of Bergstrom et al. |32| and Beyer et al. |28| the value of < v 2 > is 



found to be significantly larger for the J/ip than for the T, as intuitively expected, given 
the smaller mass of the charm quark. The near equality of the values of < v 2 > for the 
J /if) and the T and the very large value found for the T in References |52J and may 
be a consequence of an extreme choice of the parameters of the potential in the case of 
theT. 



In the present paper, more weight is given to the more recent results of Chiang et al, 
Chao et al and Schuler which are roughly consistent with each other. In the following, the 
values taken are: < v 2 >= 0.28 for the J/ip and < v 2 >= 0.09 for the T, which lie near 
the middle of the range of values obtained by these last three authors. As it will be seen 
that the effects of relativistic corrections on the shape of the fitted photon spectra are, 
after the inclusion of gluon mass effects, small (as already conjectured in Refs [15|, 0) the 



conclusions of the present work are not sensitive to the precise values assumed for < v >. 
The relativistically corrected inclusive photon spectra for the J /if) and the T calculated 
using Eqn(2.8) are shown, in comparison with the lowest order (LO) QCD prediction |33 



in Figs. la and lb respectively. Also shown in Fig. la is the curve given by truncating the 
correction to the decay rate above 0(v 2 ). In this case the spectrum is set to zero if the 
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prediction is negative. It may be remarked that Eqn(2.8) shows singular behaviour as 
z — ► 00. However, this does not affect any of the fits presented below, as no experimental 
measurements exist for z < 0.2. 



3 Higher Order QCD Corrections 



2.5 




Figure 2: QCD predictions for the inclusive photon spectrum in radiative T decays. 
Dot-dashed: LO WM, dotted: RDF ffifij , dashed: Photiadis / fT/f/ and solid curve: 
Kramer(NLO) /(X|/. 



To lowest order (LO) in perturbative QCD, the inclusive photon spectrum is described 
by the process V — > 7gg. Assuming massless gluons and neglecting relativistic corrections, 



the shape of the photon spectrum is the same as in orthopositronium decay |33 



1 dT 

r dz 



(tt 2 - 9) 



4(1 - z) ln(l 



4(l-z) 2 ln(l-^) 2z{2 



[2-z] 



z) 2 ; ( I 



(2-z) 



(3.1) 



The first estimation of higher order QCD corrections to the spectrum was made by 
R.D. Field (RDF) |j34| . These QCD effects were calculated using a parton shower Monte 
Carlo program in which the process g — > gg was iterated. The invariant mass of the cas- 
cading virtual gluons was cut off at the scale /i c = 0.45 GeV and a value A = 0.2 GeV 
was used for the QCD scale parameter in the parton shower. For T decays the average, 



The corresponding peaks near z = 0.0 are suppressed in Fig.(l) 
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perturbatively generated, 'effective gluon mass', i.e. the mass of the virtual gluon initi- 
ating the parton cascade, was 1.6 GeV. Because of the low value of the cut-off scale, the 
shapes of the photon spectra for J/ip and T decays were predicted to be similar. In both 
cases the average value of z was found to be 0.57, and even in the case of the decay of a 
hypothetical state with a mass of 60 GeV, the average z increased only to 0.59. The RDF 
spectrum for T decays is shown in Fig. 2 as the dotted line. The parton cascade used by 
RDF does not take into account QCD coherence effects in gluon radiation usually 
implemented in parton shower Monte Carlo event generators by an 'angular ordering' 



ansatz |5Bfl . The effect of this coherence, which is the QCD analogue of the 'Chudakov 
Effect' |37 in QED, is to suppress corrections due to real gluon radiation for kinematical 
configurations that yield photons close to the kinematical end point. In this case the 
two primary recoiling gluons are almost collinear, forming an effectively colourless cur- 
rent from which the radiation of large angle secondary gluons is strongly suppressed by 
destructive interference |38|]. The radiation of almost collinear gluons is not suppressed, 
but such radiation will hardly modify the shape of the LO spectrum. It is thus to be 
expected that the neglect of QCD coherence in the parton shower used to calculate the 
RDF spectrum will result in a too strong suppression of the spectrum in the end point 
region. The comparison, shown below, with a complete NLO perturbative QCD calcula- 
tion, where such coherence effects are taken into account, indicates that this is indeed the 
case. The RDF spectrum gives a good description of four out of five of the experimental 
measurements of the T spectrum (see Section 6 below). It will be seen however, that, 
for the case of J/?/> decays, the spectrum is much too hard to describe the experimental 
measurements. 

The second estimate of higher order QCD corrections to the inclusive photon spectrum 
in T decays was made by Photiadis |jT7|| . This calculation, expected to be valid near to 
z = 1, neglected completely real gluon radiation, which, as discussed above, is strongly 
suppressed in this region, but resummed to all orders in a s the leading logarithmic terms 
of the form ln(l — z) resulting from the exchange of virtual gluons and quarks between 
the two recoiling gluons of the LO diagram. As shown in Fig. 2 (the dashed curve), these 
effects give only a modest suppression of the LO spectrum near z — 1. 

The most recent result on higher order QCD corrections to the T spectrum is the 
complete NLO calculation of Kramer [ITJJ] that is also shown in Fig. 2 as the solid curve. 



It can be seen that strong suppression occurs only very near to z = 1, and is much less 
marked than in the case of the RDF spectrum. 

Since the Photiadis calculation does not include the effects of real gluon radiation, 
it can be argued that the corrections calculated by both Kramer and Photiadis should 
be applied. This will double count virtual corrections of the type shown in Fig. 2c of 
Reference |14| , but should give a better description, particularly away from the end point 



region, than using only the Photiadis correction. 

To date, no calculations of the inclusive photon spectrum in heavy quarkonia decays 
have been made taking into account, at the same time, higher order QCD corrections, 
relativistic corrections and genuine gluon mass effects []. Indeed, since the pioneering 

4 That is, including a fixed gluon mass in the calculation of both the invariant amplitude and the phase 
space, and taking into account the longitudinal polarisation states of the gluons. This is to be contrasted 
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paper of PP [12 1, gluon mass effects have been completely calculated only for the LO 



processes: V — ► ggg and V — > 7gg |39j (see Section 4 below). 

In the present analysis, the higher order QCD calculations of Photiadis and Kramer, 
made for T decays are also used, unmodified, for J/ip decays. In fact it will be seen that 
the observed end point suppression of the photon spectrum of the J/ip is so large, as 
compared to the predicted effect of both relativistic and higher order QCD corrections, 
that these play only a minor role. Indeed, the value of the effective gluon mass, m g , 
needed to describe the experimental spectrum is little affected by the inclusion of these 
corrections. The ansatz used to apply the higher order (HO) QCD corrections is to 
multiply the relativistically corrected photon spectrum given by Eqn(2.8) by the QCD 
correction factor: 

dT HO ,dT LO 

QCD ~^tz~'^tz~' { ] 

In the case of fits with m g ^ 0, phase space limitations are taken into account by the 
replacement: z — > z/zmax m Eqn(3.2) where 

zmax = 1 - -j^2~- (3.3) 

In view of the large value found for the ratio m g /Mv for the J/ip, it is to be expected 
that phase space suppression effects will be even more important for the HO corrections 
than for the LO process. This will reduce even further the effect of such corrections on 
the fitted value of m„. 



4 Gluon Mass Effects 



The possibility of gluon mass effects in the decay J/ip —> 7X was first considered by 



PP [12|. They noted that the very strong suppression of the end point region of the photon 



spectrum measured by the Mark II collaboration [13 could be explained by introducing a 
gluon mass of about 0.8 GeV. The comparison of the PP prediction with the experimental 
data did not, however, take into account experimental resolution effects which are very 
large in this case. Also relativistic and higher order QCD corrections were not included. 
The aim of the present paper is perform a similar comparison to that made by PP, but 
using all available experimental data on both J/tp and T radiative decays, as well as 
including experimental resolution effects, relativistic corrections and higher order QCD 
corrections. As the latter two effects also suppress the end point region of the photon 
spectrum, only a complete quantitative analysis, including all relevant effects, can show 
if the introduction of a non-vanishing effective gluon mass is required to describe the 
experimental data. 



with the parton shower model of RDF where an effective gluon mass (actually a time-like gluon virtuality) 
is perturbatively generated from massless gluons. The distinction between 'genuine' and 'effective' gluon 
masses is discussed further in Section 4. 
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Introducing an 'effective gluon mass' in the calculation of the Born diagram has two 
effects: 

(i) Restriction of the available phase space, i.e. modification of the boundary of the 
Dalitz plot. 

(ii) Contributions to the amplitude from longitudinal gluon polarisation states. 

As will be seen, the effect (i) is, by far, the most important. 

The gluon mass is 'effective' because it is defined only at the level of the Born diagram. 
When such a prediction is fitted to the data, which includes QCD corrections to all orders, 
the value of m g is expected to be different from the value obtained if the gluon mass were 
correctly included also in higher order diagrams in the prediction. In fact, a 'genuine' 
gluon mass M g might be operationally defined as the effective mass corresponding to a 
hypothetical all-orders pQCD calculation with massive gluons. If phase space limitations 
are very important, as in the case of the J/V', the tree level 'effective' value is expected 
to be lower than the 'genuine' value that would be found in a fit that properly includes 
gluon mass effects at all orders. If, on the other hand, the 'genuine' gluon mass is small 
compared to that of the decaying state, the effects of the non-vanishing gluon mass will 
be limited to a small region near the boundaries of phase space. In this case the tree 
level 'effective' mass is mainly generated perturbatively by splitting into gluon and quark 
pairs (as in the RDF model) and is expected to be much larger than the 'genuine' value. 
However, the 'genuine' mass M g found by comparing the prediction of the all-orders pQCD 
calculation to the data is expected to be independent of the mass of the decaying state. 
The above argument also shows that, for some mass of the decaying state, the tree level 
'effective' mass and the 'genuine' gluon mass should be equal. It may be conjectured that 
this is almost the case for the T. 

The correction curves for the processes J ftp — > ggg and rj c —>■ gg calculated by PP 
and shown in their Fig.l took into account both the effects (i) and (ii), but no explicit 
formulae were given. In the present paper essential use is made of formulae including 
both effects (i) and (ii) obtained by Liu and Wetzel f3"9|| . For the decay V — > 7gg, the 
fully differential spectrum is: 

1 dT 1 1 . 8 , 25 w , 2 



T dzdx x dx 2 (tt 2 - 9) z 2 (x' 2 ) 2 (x' 3 ) 2 3 2 

+ 32,2(1 - 2r])(l - |)r? 2 + z 2 {l - 2rj)(l - Gq - Qrf) 

- 2(z 3 + (x 2 ) 3 + (4) 3 )(1 - f V - 2 V 2 ) + z\\ + |) 

+ ((4) 2 + 04) 2 )(1 - 8^ + 22r? + 8r/ 3 ) - (x' 2 ) 2 {x' 3 ) 2 r, 

+ (^ 4 + (x 2 ) 4 + (4) 4 )(l + |)] (4.1) 



2E 7 /My, T] = (m g /Mv) 2 , Xi = 2Ei/My, x\ = Xi — 2r/, Ei = gluon energy 
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Experiment 


Decay process 


K = <W# 7 ( GeV ) 


/b 7 at £ 7 = M v /2 (%) 


Mark II 


J/V> -> 7X 


0.12/jRy 


9.6 


i)' -> 7T + 7T-(J/^ -> 7 X) 


, 1 

0.12/ y/Ey 


9.6 


J/V> -" (7 e+e")X 


0.022E^ 2b 


3.4 


CUSB 


T — > 7X 


0.039/£°- 25 


3.2 


ARGUS 


T — > 7X 


^0.0052 + 0.0042/£ 7 


7.8 


Crystal Ball 


T — > 7 A" 


0.027/E^ 5 


1.8 


CLE02 


T — > 7 A 


0.0035/£0- 75 + 0.019 - 0.001£ 7 


1.5 



Table 2: Photon energy resolutions of the different heavy quarkonia radiative decay exper- 
iments. 



Here To is the radiative width uncorrected for gluon mass effects. The allowed phase 
space region is defined by the conditions^: 

(4.2) 
(4.3) 
(4.4) 



2 = z + x 2 + x 3 
< z < 1 - 4r? 

c™ in < x 2 < x™ ax 



max 
J/ 2 



mm 

x 2 




(4.5) 
(4.6) 



In the gluon mass dependent fits to be presented below, the functions fo(z) and fi(z) 
in the KM formula (2.8) are replaced by the functions: 



h(z,m g ) 



dx 2 f (z,x 2 ,m g ) 
dx 2 fi{z,x 2 ) 



(4.7) 
(4.8) 



Where fo{z,x 2 ,m g ) is derived from Eqn.(4.1) and f\(z,x 2 ) from Eqn.(3.5) of KM. Thus 
phase space effects are taken properly into account in both f and fx, whereas the effect 
of longitudinal polarisation states are included only in f Q . As will be shown below, the 
latter effect is much smaller than the former, so that the effect on the fit results of 
the uncalculated contribution of gluon longitudinal polarisation states on the relativistic 
correction coefficient f\ is expected to completely negligible. 



5 The Decay J/tp -> 7X 



To date, the inclusive photon spectrum in J/i/j decays has been measured by only one 
experiment, the Mark II collaboration ||13|| . Actually, in Reference ||13|| , three independent 



3 Note that r/ is defined differently here than in Reference (l( 
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Figure 3: Fits, assuming m g = 0, to inclusive photon spectra in J/ip decays. Mark II (top 
left), Mark II cascade (top right), Mark II conversion (bottom). Dash-dotted line: LO 
QCD prediction dashed line: Relativistic Correction (RC) included J^HjJ , solid line: 
RC and Photiadis [TjJ HO QCD correction, fine dotted line: RC and Kramer t l^J HO 
QCD correction and dotted line: RC, Photiadis and Kramer corrections 
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Figure 4: Fits, for the effective gluon mass m g , to inclusive photon spectra in J /if) decays. 
Mark II (top left), Mark II cascade (top right), Mark II conversion (bottom). Fit curves 
are defined as in Fig. 3. 
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Channel C 


z 


BRxlO 4 


Relative BR 


Resonances 


7T° 


0.998 


0.39 ± 0.13 


0.009 


77 

1 


0.969 


8.6 ± 0.8 


0.2 


n' 


0.904 


43.1 ± 3.0 


1.0 


77(1440) -> op 


0.784 


17.0 ± 4.0 


0.39 


77(1440) -> P7 


0.784 


0.64 ± 0.14 


0.015 


7/(1440) -> XFtt 


0.784 


9.1 ± 1.8 


0.21 


7/(1440) -> 7/7r + 7r" 


0.784 


3.0 ± 0.5 


0.07 


77(1760) -> pp 


0.677 


1.3 ± 0.9 


0.030 


772(1870) -> 7T7T 


0.63 


6.2 ± 2.4 


0.144 


7/(2225) 


0.484 


2.9 ± 0.6 


0.067 


/ 2 (1270) 


0.832 


13.8 ± 1.4 


0.32 


/i(1285) 


0.828 


6.1 ± 0.9 


0.14 


/i(1420) -> AT/frr 


0.790 


8.3 ± 1.5 


0.19 


/o(1500) 


0.765 


5.7 ± 0.8 


0.13 


/l(1510) ^T/TT + TT- 


0.762 


4.5 ± 1.2 


0.10 


#(1525) 


0.758 


4.7 ± 0.6 


0.11 


/o(1710) - 


0.695 


8.5 ± 1.0 


0.20 


/ 2 (1950) -> K*^* 


0.603 


7.0 ± 2.2 


0.16 


/ 4 (2050) 


0.561 


27.0 ± 7.0 


0.63 


/j(2220) -> 7T7T 


0.486 


0.8 ± 0.4 


0.019 


/j(2220) -> 


0.486 


0.8 ± 0.3 


0.019 


/j(2220)^pp 


0.486 


0.15 ± 0.08 


0.003 


Total BR (Res.) 


180 ± 10 


4.17 


Exclusive states 


7T + 7T 7r°7r° 




83.0 ± 31.0 


1.93 


7T + 7T _ 7r + 7r _ 




28.0 ± 5.0 


0.65 






21.0 ± 6.0 


0.49 


777T7T 




61.0 ± 10.0 


1.42 


PP 




45.0 ± 8.0 


1.05 


UJijJ 




15.9 ± 3.3 


0.37 


00 




4.0 ± 1.2 


0.093 






40.0 ± 13.0 


0.92 






3.8 ± 1.0 


0.09 


Total BR (Excl 


•) 


302 ± 37 


7.0 


Total BR (Res.+Excl.) 


482 ± 38 


11.2 



Table 3: Composition of the hadronic final state in the 3 /if) radiative decays: J /if) 
7C# 
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7]' fraction 


7] fraction 


m g 


v 2 . 
Amm 


CL 


Mark II (25 dof) 


0.0 


0.0 


0.610 ± 0.015 


51.7 


2xl0~ 3 


0.0741°,;™? 


0.015 


0.720 ± 0.012 


17.5 


0.86 


Mark II cascade (25 dof" 




0.0781°;°°? 


0.016 


n 799+0.U19 


18.3 


0.79 


Mark II conversion (22 dof) 


u - u ' "J-0.010 


0.015 


0.597 ± 0.019 


20.0 


0.58 



Table 4: Results of fits to determine the fractions of J/ip — *> 777', 777 in J/ip radiative 
decays 



Fitted Model 


V 2 . 

Amm 


CL 


Mark II (27 dof) 


LO 


666 


< IO- 30 


Rel. Corr n 


336 


< IO- 30 


Rel. Corr n , QCD(P) 


279 


< io- 3U 


Rel. Corr n , QCD(K) 


253 


< IO- 30 


Rel. Corr n , QCD(PxK) 


224 


< IO- 30 


Mark II cascac 


e (26 c 


lof) 


LO 


348 


< IO- 30 


Rel. Corr" 


183 


1.5xl0" 25 


Rel. Corr™, QCD(P) 


153 


5.9xl0- 20 


Rel. Corr n , QCD(K) 


136 


7.3xl0" 17 


Rel. Corr", QCD(PxK) 


121 


3.3xl0" 14 


Mark II conversion (24 dof) 


LO 


198 


< 2.5xl0~ 29 


Rel. Corr n 


80.9 


4.4xl0" 25 


Rel. Corr n , QCD(P) 


66.3 


7.8xl0" 6 


Rel. Corr n , QCD(K) 


50.8 


l.lxlO" 3 


Rel. Corr n , QCD(PxK) 


45.0 


5.8xl0" 3 



Table 5: Results of fits with m g = to J/ip inclusive photon spectra. LO denotes the lowest 
order Q CD prediction (Eqn. (3.1)). Rel. Corf 1 includes relativistic corrections calculated 
according to Eqn. (2.8). QCD(P) and QCD(K) denote, respectively the Photiadis j[F% and 
Kramer HO QCD corrections. QCD(PxK) means that both corrections are applied. 
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Fitted Model 


m g (GeV) 


V 2 . 

Amm 


CL 


Mark II (26 dof) 


LW 


0.734 ± 0.010 


24.6 


0.54 


LW, Rel. Corr n 


n 74 n+0.0Uy 


22.6 


0.66 


LW, Rel. Corr n 








QCD(P) 


n 791 +0.010 

u - ' zi -0.009 


17.5 


0.89 


LW, Rel. Corr n 








QCD(K) 


0.6651H1! 


19.9 


0.80 


LW, Rel. Corr n 








QCD(PxK) 


0.6531S 


16.9 


0.91 


Mark 1 


I cascade (25 dof) 


LW 


D 7Q7+0.015 
u - ' ° ' -0.018 


24.4 


0.50 


LW, Rel. Corr n 


0.740± 0.017 


21.9 


0.64 


LW, Rel. Corr n 








QCD(P) 


0.719± 0.019 


18.4 


0.82 


LW, Rel. Corr n 








QCD(K) 


0.667± 0.021 


22.7 


0.60 


LW, Rel. Corr n 








QCD(PxK) 


0.655± 0.021 


20.0 


0.75 


Mark II 


conversion (23 dof) 




LW 


0.623™ 


31.6 


0.11 


LW, Rel. Corr n 


0.607+H^ 


23.9 


0.41 


LW, Rel. Corr n 








QCD(P) 


0.598l» 


19.9 


0.65 


LW, Rel. Corr n 








QCD(K) 


n C07+0.025 

U.O<3 ( _o.030 


22.5 


0.49 


LW, Rel. Corr n 








QCD(PxK) 


5?fi +0 - 027 

U.OZO_g 29 


20.3 


0.62 



Table 6: Results of fits with variable m g to J/i/j inclusive photon spectra. LW denotes 
gluon mass corrections calculated according to the calculations of Liu and Wetzel JjS^jJ . 
Relativistic and HO QCD corrections are defined as in Table 5. 
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measurements of the spectrum are given. The first (referred to simply as 'Mark IF) uses 
the process: 

e + e~ — > J /if) — ► 7X 

where the photons are detected in the electromagnetic calorimeter of the detector. The 
second, sample, 'Mark Il(cascade)', uses the process: 

e + e~ — > if)' — > J/lpTTTC — > 7X7T7T 

In this case the acceptance and resolution are similar to the 'Mark IF measurement. 
The third data sample, 'Mark II (conversion)', uses J/?/> radiative decay events where the 
photon converts into an e + e~ pair in the beam pipe or the inner flange of the tracking 
chamber. Measurement of the momenta of the e + ,e~ in the chamber yields a sample 
with reduced statistics but much improved photon energy resolution The photon energy 
resolutions for these three event samples are given in Table 2. 

As well as these inclusive measurements, many exclusive measurements have been 
made where a single resonant state or an exclusive multihadron final state is produced 
in association with a hard photon ||. These measurements are summarised in Table 3 
where the values of z for each exclusive channel with a single particle recoiling against 
the photon are given. As can be seen from Table 3 ( see also Fig. 59 of Reference |[40|| ) the 
most striking feature of the photon spectrum near the end point is the strong exclusive 
production of 77 and rf mesons: 

J/<0->77?, 77/ 

These two channels alone account for 13% of the total branching ratio for z > 0.6 and 
completely dominate the end point region z > 0.85. Because of the large contribution of 
these two resonances, it was not possible to obtain acceptable fits to the spectra using 
the function of Eqn.(2.8) even when HO QCD corrections and gluon mass effects were 
included. The procedure adopted was then to fix the ratio r(J/^> — > ■yr]')/r(J /ip — > 777) 
to the measured value 5.0 ||, and perform fits treating the ratio: 

R v > = r(J/V — > 77/)/T , (J/'0 — > 7 continuum) 

as a free parameter, which include the 77' and 77 contributions at the appropriate z values of 
0.904 and 0.969 respectively. Here '7 continuum' refers to Eqn(2.8), including also gluon 
mass effects and the Photiadis HO QCD correction. The other two parameters in the fit 
are an overall normalisation constant and the effective gluon mass m g . The method used 
to fold in the experimental resolution function is described in the Appendix. As shown 
in Table 4, attempting to fit the spectrum without explicitly introducing the 77' and 77 
contributions (i.e. with = 0), leads to an unacceptably low confidence level (CL) of 3 
x 10~ 3 . However, including their contributions, good fits are obtained for all three spectra 
with consistent values of R v i . Their weighted average is: 

R v , = 0.0754 ± 0.0070 

The Mark II and Mark Il(cascade) spectra yield consistent values of m g around 720 MeV, 
but the Mark II (conversion) spectrum gives a significantly lower (5.8 a) value of 597 MeV, 
indicating some systematic difference in the latter measurement. The weighted average 
value of R^i corresponds to : 

BR{J/ip -> 7 + continuum + 77 + 77') = 0.063 ± 0.0070 
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which may be compared to the summed branching ratio of all channels reported in Table 
3 of 0.0482 ± 0.0038. So it is estimated that 76% of all J/ip radiative decays to light 
hadrons are contributed by the exclusive channels listed in Table 3. 

The results of fits to the three inclusive photon spectra, including the exclusive rj' 
and i] contributions, estimated using the weighted average value of obtained above, 
without including gluon mass effects, are presented in Table 5 and Fig. 3. No acceptable 
fits are obtained, even after the inclusion of of both relativistic and HO QCD corrections. 
The best confidence levels obtained for the fits to the Mark II, Mark II (cascade) and Mark 
II (conversion) spectra are: < 10~ 30 , 3.3 x 10~ 14 and 5.8 x 10 -3 respectively. Although, as 
expected from Fig. 1, the shape of the spectrum is drastically modified by the relativistic 
correction, the change in shape, though qualitatively in the right direction, is, by far, not 
enough to explain the observed spectrum shape. As can be seen in Fig. 3a, the estimated 
effects of HO QCD corrections are even smaller than those of the relativistic correction. 

A similar series of fits, but including gluon mass effects, is presented in Table 6 and 
Fig. 4. Good fits are obtained in all cases, and it can be seen that the inclusion of rel- 
ativistic and HO QCD corrections have only a minor effect on the fitted values of m g . 
For example, in the case of the largest statistics (Mark II) data sample, introducing the 
relativistic correction increases the fitted value of m g of 734 MeV by only 6 MeV. Further 
applying either the Photiadis or Kramer HO QCD corrections give further shifts of -19 
MeV or -75 MeV respectively in the fitted m g value. In fact because of the relatively large 
value of m g as compared to Mj/^, it is to be expected, as mentioned previously, that HO 
QCD corrections will be much reduced by by phase space and propagator suppression 
effects. Choosing then the fit with relativistic and the smallest (Photiadis) HO QCD 
correction as best estimate yields the result: 

m g = 0.7211°;™ t° ° ll GeV (Jyty) 

where the first error is statistical (from the fit to the Mark II spectrum) and the second is 
systematic, conservatively estimated from the full spread of the different fit results to the 
Mark II data given in Table 6. In view of the size of this systematic error in the Mark II 
value, and the large observed systematic shift in m g obtained with the Mark II (conversion) 
spectrum, no significant improvement in the knowledge of m g is expected by combining 
the results of the fits to the three spectra. The less precise values provided by the Mark 
Il(cascade) and Mark Il(conversion) data should then be considered as consistency checks. 

It may be noted that the exclusive rj' signal is clearly seen in the Mark II (conversion) 
spectrum shown in Figs. 3c and 4c. The shape of the observed peak is well described 
in Fig. 4c by the fit including the relativistic and Photiadis HO QCD corrections, and 
the experimental resolution function given in Table 2. This agreement gives an important 
cross check on the method used here (see the Appendix) to fold in the experimental 
resolution effects. 

In order to study the relative importance of phase space effects and longitudinal gluon 
contributions in the gluon mass fits, the fit to the Mark II data, including both the 
relativistic and the Photiadis HO QCD corrections, is repeated setting m g = in the 
function f (z,X2,m g ) of Eqn(4.7), thus removing the longitudinal contributions. A very 
good fit is still obtained (confidence level 0.94) with m g = 0.682 + 0.010-0.013 GeV. This 
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Fitted Model 


A mm 


CL 


CUSB (11 dof) 


LO 


30.2 


< 1.5xl0~ 3 


Rel. Corr™ 


10.9 


0.45 


Rel. Corr™, QCD(P) 


14.7 


0.2 


Rel. Corr™, QCD(K) 


13.0 


0.29 


Rel. Corr n , QCD(PxK) 


21.3 


0.03 


ARGUS (19 dof) 


LO 


95.3 


3.8xl0" 12 


Rel. Corr n 


62.6 


1.5xl0" 6 


Rel. Corr™, QCD(P) 


44.7 


7.5xl0" 4 


Rel. Corr™, QCD(K) 


49.8 


1.4xl0" 4 


Rel. Corr n , QCD(PxK) 


37.6 


6.7xl0~ 3 


Crystal Ball 


[14 dof) 


LO 


51.9 


2.9 10~ 6 


Rel. Corr n 


31.6 


4.6xl0~ 3 


Rel. Corr™, QCD(P) 


22.4 


0.071 


Rel. Corr™, QCD(K) 


24.6 


0.039 


Rel. Corr n , QCD(PxK) 


18.9 


0.17 


CLE02 (22 dof) 


LO 


90.7 


2.6xl0" 10 


Rel. Corr n 


56.9 


6.3xl0" 5 


Rel. Corr™, QCD(P) 


32.1 


0.076 


Rel. Corr™, QCD(K) 


48.8 


8.5xl0" 4 


Rel. Corr n , QCD(PxK) 


29.5 


0.131 



Table 7: Results of fits with m g = to T inclusive photon spectra. See Table 5 for the 
definitions of the different fits. 



is only 5.4% lower than the value m g = 0.721 + 0.010 — 0.009 GeV obtained including the 
longitudinal contributions. Thus the gluon mass correction is dominated by phase space 
effects. 



6 The Decay T —> 7X 



Five different experiments have measured the inclusive photon spectrum in T decays: 
CUSB 



H, CLEO H, ARGUS |3 



p. The CLEO 



Crystal Ball @ and CLE02 
measurement, which, like Crystal Ball, ARGUS and CLE02, but unlike CUSB, is in good 
agreement with the RDF prediction, is not analysed in the present paper as no efficiency 
corrected spectrum was provided. 

Unlike for the case of the J/ip, no positive evidence has been found for the exclusive 
production of single resonances in the radiative decays of the T ||. The efficiency cor- 
rected inclusive photon spectra measured by CUSB, ARGUS, Crystal Ball, and CLE02 
have therefore been directly fitted to the relativistically corrected spectrum Eqn(2.8), 
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Figure 5: Fits, assuming m g = 0, to inclusive photon spectra in T decays: CUSB (top 
left), ARGUS (top right), Crystal Ball (bottom left), CLE02 (bottom right). Fit curves 
are defined as in Fig. 3. 
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Figure 6: Fits, for the effective gluon mass m g , to inclusive photon spectra in T decays: 
CUSB (top left), ARGUS (top right), Crystal Ball (bottom left), CLE02 (bottom right). 
Fit curves are defined as in Fig. 3. 
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Fitted Model 


m g (GeV) 


Amtn 


CL 


CUSB (10 dof) 


LW 


0.66 ± 0.08 


5.2 


0.88 


LW, Rel. Corr n 


0.64 ± 0.09 


6.1 


0.81 


LW, Rel. Corr n 








QCD(P) 


0.54 ± 0.12 


13.1 


0.22 


LW, Rel. Corr n 








QCD(K) 


o.i6i^ 6 7 


12.6 


0.25 


LW, Rel. Corr n 








QCD(PxK) 


o.i5±8:ii 


21.0 


0.021 


ARGUS (18 dof) 


LW 


1 on+U.U8 
1 - oy -0.10 


27.3 


0.074 


LW, Rel. Corr n 


1.39± 0.10 


28.8 


0.051 


LW, Rel. Corr n 








QCD(P) 


-i 97 +0.11 
L - Zl -0.12 


23.3 


0.18 


LW, Rel. Corr n 








QCD(K) 




34.4 


0.011 


LW, Rel. Corr n 








QCD(PxK) 


i.06l8il 


28.8 


0.05 


Crystal Ball (13 dof) 


LW 


i o i+O.lO 
i - zi -0.09 


21.5 


0.064 


LW, Rel. Corr n 


i 91+0.10 
i - z± - 0.11 


22.2 


0.052 


LW, Rel. Corr n 








QCD(P) 


i-i4±g : g 


18.0 


0.16 


LW, Rel. Corr n 








QCD(K) 


0.90i8i 


23.5 


0.036 


LW, Rel. Corr n 








QCD(PxK) 


No y 2 • 






CLE02 (21 dof) 


LW 


1.27 ± 0.07 


29.5 


0.103 


LW, Rel. Corr n 


i 9 c;+0.U8 


30.4 


0.08 


LW, Rel. Corr n 








QCD(P) 


i i c;+0.08 
i - id -0.09 


16.9 


0.72 


LW, Rel. Corr n 








QCD(K) 


I nq + 0.09 

-'-•U'J-0.12 


34.5 


0.032 


LW, Rel. Corr n 








QCD(PxK) 


o.9o±8:ii 


22.5 


0.37 



Table 8: Results of fits with variable m g to T inclusive photon spectra. The descriptions 
of the different fitted models are the same as in Table 6. 
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possibly also including HO QCD corrections as discussed in Section 3 and gluon mass 
effects as described in Section 4. Experimental resolution effects are included in the same 
way as for the fits to the J/ip decays described above. The photon energy resolutions 
of the different experiments are given in Table 2. Results of fits assuming a vanishing 
effective gluon mass are presented in Table 7 and Fig. 5, while fits to m g and an overall 
normalisation constant yield the results shown in Table 8 and Fig. 6. 

As can be seen in Figs. 5 and 6, the CUSB spectrum differs markedly in shape from 
those of the later ARGUS, Crystal Ball, and CLE02 experiments. The suppression of 
the end-point region, relative to the LO QCD prediction (the dash-dotted curves), is 
much reduced. The results of the fits to the CUSB data should be treated with caution 
as, unlike the other experiments, the published errors on the photon spectrum are are 
purely statistical. A relatively large systematic error is expected (as found in the other 
experiments) from the 7r°/7 separation procedure, especially at lower values of z. 

Considering first the fits with m g = in Table 7, it can be seen that the LO spectrum 
is ruled out, individually, by all four experiments. Inclusion of the relativistic correction 
with < v 2 >= 0.09 gives a good description of the CUSB data but is ruled out with 
a confidence level of less than 0.5% by each of the other experiments. The best overall 
description is given by combining the relativistic and Photiadis HO QCD corrections. 
Low, but acceptable, confidence levels of 0.2, 0.071 and 0.076 are found for the fits to 
CUSB, Crystal Ball and CLE02 data, respectively. Only the ARGUS spectrum (CL = 
7.5xl0~ 4 ) is inconsistent with this hypothesis. However, combining the fits to all four 
experiments gives x 2 /dof = 114/66 (CL=2.3x 10~ 4 ). 

For ARGUS, Crystal Ball and CLE02 the best fits are given by including both the 
Photiadis and Kramer HO QCD corrections. The combined fit gives, however, x 2 /dof = 
86/66 (CL=4.7x 10 -3 ). Thus no consistent overall description of the data is found for 
m g = 0. 

When the effective gluon mass m g is included as a fit parameter it can be seen (Table 8) 
that fits with confidence levels > 1% are found for all fit hypotheses and all experimental 
spectra. However, the ARGUS, Crystal Ball and CLE02 data give consistent values of m g 
in the range 0.9-1.4 GeV, whereas significantly lower values 0.15-0.66 GeV are found in the 
fits to the CUSB spectrum. Because of this discrepancy and the neglect of (potentially 
large) systematic errors in the latter experiment, only the results from the three most 
recent experiments are used to obtain the average value of m g quoted below. For these 
experiments the best overall fit is given (as in the case m g = 0) by including both the 
relativistic and the Photiadis HO QCD correction. This yields, for the weighted average 
value of the effective gluon mass: 

m g = 1.18 ±0.06 GeV 

where the error quoted is derived from fit errors of the different experiments. Performing 
a fit to the ARGUS, Crystal Ball and CLE02 spectra, with m g fixed at the above value, 
and varying only the normalisation constants of the fitted curves gives a good overall fit 
with \ 2 /dof = 59.0/55 (CL=0.33). Making the same type of fit, but including CUSB 
leads to: x 2 /dof = 121.0/66 (CL=4.3xlO~ 5 ). The published CUSB data is therefore 
clearly inconsistent with value of m g favoured by ARGUS, Crystal Ball and CLE02. 
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This apparent inconsistency of the CUSB measurement, is, however, very sensitive to 
the error assignement of the data. Increasing the quoted (purely statistical) errors by a 
constant factor of 1.5 to account for systematic effects, modifies the last fit result quoted 
above to: x 2 /dof = 89.0/66 (CL=0.031). The CUSB measurement is now marginally 
consistent with the average of the three other experiments. 

The theoretical systematic error on m g is estimated in the same way as for J ftp decays. 
This gives: 

m g = 1.18 ± 0.06lg;^ GeV (T) 

where the first error is experimental and the second is a conservatively estimated theo- 
retical error that that includes the full range of relativistic and HO QCD corrections in 
the fits to the CLE02 data in Table 8. 

To investigate the importance of the effects of longitudinal gluon polarisation states 
for the case of the T, a fit is made to the CLE02 data including relativistic and the 
Photiadis HO QCD corrections, with m g = in the function fo of Eqn(4.7). The fit gives 
m g = 1.10 + 0.08 - 0.09 GeV with x 2 /dof = 13.1/21 (CL=0.91) to be compared with 



ni 



•A 



1.15 + 0.08 - 0.09 GeV with x 2 /dof = 16.9/21 (CL=0.70) when longitudinal gluon 
contributions are included in f . As in the case of the J/ip, longitudinal gluon states give 
only a small effect; they increase the fitted value of m g by only 4.5%. 



Finally in this Section a few remarks on the CLEO measurement ji2], which is not 
included in the present analysis. Relativistic corrections were not taken into account in the 
theoretical predictions, which were suitably modified to account for detector acceptance 
and resolution effects before comparison with the (uncorrected) experimental data. Fitting 
the LO QCD, Photiadis and RDF spectra to the data yields \ 2 °f 14.2, 10.0, and 8.1, 



respectively, for eleven degrees of freedom p2] . The corresponding respective confidence 



levels are 0.22, 0.46 and 0.70. Although the RDF model is slightly favoured, all the 
fits have acceptable confidence levels and no distinction between the different theoretical 
models is possible from this measurement. 
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7 Determinations of a s (rriQ) 




QGeV 

Figure 7: Comparison of a s (m c ) and a s (mb) values obtained from measurements of R'jm 
and R' r compared to the world average value (dashed curve with ±1<t limits indicated by 
the solid curves) from Reference /Q/. The crosses with error bars show values obtained 
assuming m g = and i^/tuq = 1.0. The effect of varying the renormalisation scale in the 
range 0.6 < /V m Q < 2.0 is indicated by the vertical boxes. The square points with error 
bars are obtained by applying full gluon mass corrections with m g = 0.721 GeV for the 
J/ip and m g = 1.18 GeV for the T with n/rriQ = 1.0. NLO QCD corrections are applied 
in all cases. 
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Branching ratio 


Experimental value 


Reference 


r(J/0 — > hadrons) /rj/,/, 


0.632 ± 0.022 




r(J/0 — > 7 + hadrons) /Tj/,/, 


0.0624 ± 0.0067 


This paper 




10.13 ± 1.14 


Ratio of above 




33.33 ± 2.44 


ARGUS g3 






37.04 ± 6.17 


Crystal Ball [ 


g 




36.36 ± 2.11 


CLE02 45 






35.46 ± 1.51 


Weighted average 



Table 9: Experimental branching ratios used to determine a s (mQ) 





0.6 


1.0 


2.0 


a s (m c ) 
A(PDG) 
Deviation(cr) 


0.178 ± 0.011 
-0.179 ± 0.022 
-8.1 


0.184 ± 0.015 
-0.173 ± 0.024 
-7.2 


0.213 ± 0.023 
-0.144 ± 0.03 
-4.8 


a s {m b ) 
A(PDG) 
Deviation(cr) 


0.169 ± 0.005 
-0.048 ± 0.008 
-6.0 


0.169 ± 0.006 
-0.048 ± 0.009 
-5.3 


0.186 ± 0.008 
-0.031 ± 0.01 
-3.0 



Table 10: q s (i)1q) values obtained neglecting gluon mass corrections for different choices of 
the renormalisation scale /i. A(PDG) is the difference from the PDG average values: 
a s (m c ) = 0.357tHiL oi s (m b ) = 0.2171^07- 'Deviation' is A(PDG) divided by its error 

The strong coupling constant a s {mq) may be determined from the experimental mea- 
surements of the branching ratio: 

T(V -> hadrons) 
v T(V ^7 + hadrons) 1 ' ' 

where V denotes a vector heavy quarkonium ground state (J/0 or Y). Use of R' v has the 
advantage, as compared with other branching ratios sensitive to a s (for example Ry = 
T(V — + hadrons)/r(U — > £ + £~), that relativistic corrections cancel in the approximation 
used in the present paper. 

In Eqn(7.1) the process l V — > 7 + hadrons' is understood to be the heavy quark anni- 
hilation process into light hadrons, for which the lowest order QCD process is V — >• jgg. 



J/A 






/ggg 


/ggg/ /ygg 


Phase Space 


40 +auv 

u -^ u -0.02 


o.is±8:S 


u -^°-0.02 


LW 


^2 +u - Ub 

u - dz -0.01 


o.is±8:S 


u - od -0.02 


T 


Phase Space 


0.741^ 


61 +U ib 


n 00+U.U8 


LW 


0.80^ 


61 +u ib 

u - oi -0.06 


77+ u - n 

u - ' ' -0.05 



Table 11: Gluon mass correction factors. Phase Space' indicates that gluon mass effects 
are taken into account only in the kinematic limits. 'LW means that the complete calcu- 
lation of Reference fiSQ] is used The errors quoted are derived from the total uncertainties 
in the fitted values of m g . Note that f lgg = f ggg = 1 for m g = 0. 
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/V m Q 


0.6 


1.0 


2.0 


No 0(a s ) Corr. 


a s (m c ) 
A(PDG) 
Deviation(a) 


n 991 +0.002 

-0.136 ± 0.027 
-5.0 


0.298™ 
-0.059 ± 0.025 


0-467™ 
0.110 ± 0.120 
0.9 


0-525™ 
0.168 ± 0.133 


a s {m b ) 
A(PDG) 
Deviation(a) 


0.189™ 
-0.028 ± 0.009 
-3.1 


0-193™ 
-0.024 ± 0.011 
-2.2 


n 991 +U.U35 
U.ZZl_o.oi5 

0.004 ± 0.016 
0.25 


0.248™ 
0.032 ± 0.025 
1.2 



Table 12: a s (mQ) values obtained including phase space gluon mass corrections for differ- 
ent choices of the renormalisation scale fi. See Table 10 for the definitions of A(PDG) 
and 'Deviation'. The quoted errors include the effect of the uncertainties in the gluon 
mass correction factors. 



/V m Q 


0.6 


1.0 


2.0 


No 0(a s ) Corr. 


« S K) 

A(PDG) 
Deviation(cr) 


U.ZZ^_oo05 

-0.133 ± 0.019 
-7.0 


n QQO+0.014 
U.OOZ_o,042 

-0.025 ± 0.024 
-1.0 


n 17+0.088 
u - 0i ' -0.178 

0.260 ± 0.178 
1.5 


0.681;™ 
0.324 ± 0.187 
1.7 


a s {m b ) 
A(PDG) 
Deviation(cr) 


n i q 7 +U.UU5 
U-J-y '-0.016 

-0.020 ± 0.009 
-2.2 


n 9QC-+0.010 

-0.012 ± 0.012 
-1.0 


0-239^ 
0.022 ± 0.033 
0.67 


n 971 +0.021 

u - z ' l -omo 
0.054 ± 0.040 

1.4 



Table 13: a s (mQ) values obtained, including full gluon mass corrections from Refer- 
ence for different choices of the renormalisation scale /i. See Table 10 for the defini- 
tions of A(PDG) and 'Deviation'. The quoted errors include the effect of the uncertainties 
in the gluon mass correction factors. 



Experiment 


otsinBLM) = a s (1.5GeV) 


CUSB 
ARGUS 
Crystal Ball 
CLEO 


U.ZZO_ .042 

0.225 ±0.011 ±0.019 
0.25 ±0.02 ±0.04 

n 97+0.03 +0.03 
u - z '-0.02 -0.02 


Weighted Mean 


0.228 ±0.016 ±0.011 



Table 14: Published values of a s at the BLM scale determined from T radiative decays. 
For CLEO, the measured spectrum is extrapolated using the RDF prediction. The first 
error quoted in each case is statistical, the second systematic. 
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Thus the non-annihilation process V — > 7?7q, where 77q is the lowest-lying pseudoscalar 
heavy quarkonium ground state, is not included. Similarly, l V — > hadrons' is understood 
to be the direct (strong interaction) process for which the lowest order QCD process is 
V — ► ggg. The contribution of the radiative process: V —>■ 7* — > qq — >• hadrons (branching 
ratio 17% for the J/VO is, therefore, not included. 

The experimental branching ratios used here to extract a s {mq) are summarised in 
Table 9. For the J/ip the branching ratio T(3/ij} — > 7 + hadrons) /Tj/^, where Tj/^, 
denotes the total width of the J/ip, is obtained from the fits to the branching ratio R v > 
presented in Table 4. The measured exclusive branching ratio into 777' ||, is used to 
derive, from R n i, the branching fraction for '7 continuum' (see Section 5 above). The 
measured exclusive 777 and 777' fractions are then added to the '7 continuum' fraction to 
give the full branching fraction into 7 + hadrons quoted in Table 9. The errors on this 
quantity are derived from that on the weighted average value of R v ' and the experimental 
777 and 777' fractions The branching fraction T(J/ip — > hadrons) /Fj/^ is derived from 
the total hadronic width of the J/ip given in Ref by subtracting the contribution || of 
the process J/tp — > 7* — > qq — ► hadrons. 

In the case of the T measurements, the branching ratio R' r was directly measured 
by CUSB, CLEO, ARGUS, Crystal Ball and CLE02. The values of R' r obtained by the 
last three of these experiments are reported in Table 9. In all cases the extrapolation of 
the measured photon spectrum to z = was done using the RDF theoretical spectrum. 
As will be shown below, the shape of this spectrum is in good agreement with the fit 
curves obtained in the present paper, which take into account relativistic corrections and 
explicit gluon mass effects at the born level. This shape is little affected by including HO 
QCD corrections according to Ref jL7| and/or Ref [fL4}| . To extract a s (m&), the weighted 
average (also reported in Table 9) of the R' r measurements of ARGUS, Crystal Ball and 
CLE02 is used. 

Consistent results for R' T were found by CUSB: 33.4 ± 6.6 and CLEO (using the RDF 
spectrum): 39.4 ± 3.6. Since, however, the shape of the inclusive photon spectrum mea- 
sured by CUSB is inconsistent with those measured by the other four experiments, and the 
analysis of the CLEO photon spectrum could not be performed, these two measurements 
are omitted from the R' r average used here to determine a s (m&). Thus, the a s analysis 
is performed using only data from experiments for which a consistent determination of 
the effective gluon mass was possible. Including also the CUSB and CLEO measurements 
in the weighted average of R' r gives the value 35.94 ± 1.36 which differs by only 0.35<r 
from the weighted average of ARGUS, Crystal Ball and CLE02 quoted in the last row of 
Table 9. 
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0.25 



0.5 



0.75 



Figure 8: Allowed regions of the Dalitz plots for V — > 7gg (open areas) and V — > ggg 
(cross-hatched areas), corresponding to the effective gluon mass values 0.721 GeV (J/ip), 
a), and 1.18 GeV (T). b). For vanishing effective gluon mass, the full areas above the 
dotted lines are kinematically allowed 
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Taking into account NLO QCD radiative corrections [46] as well as gluon mass cor- 
rections, a s {mq) may be derived from the formulae: 



Sa.M/gg [l + ^(|A)ln£-3.74)] 
16a / 7gg [i + ^^ln^-6.68)] ' 



5a^)/ ggg [l + ^(|A»ln^-4.90)] 
4a / 7gg [i + ^M^in^-r^)] ' 

Here (3q is the one-loop QCD beta function coefficient: 

where the number of active quark flavours, n./, is taken to be 3 for the J/ip and 4 for the 
T. The values of the heavy quark masses are assumed to be m c = 1.5 GeV and m& = 4.9 
GeV. The parameter \i is an arbitary renormalisation scale, and / ggg , / 7gg are correction 
factors taking into account gluon mass effects. As previously pointed out, relativistic 
corrections cancel in the ratio R' v . For any given value of the renormalisation scale fi, 
Eqns.(7.2) or (7.3) are solved for a s (/i). The corresponding value of a s {mo) is then found 
by use of the one-loop QCD evolution formula: 

1 (7.4) 



a s (m Q ) a s (n) 2tt m Q 



Setting / ggg = / 7gg = 1 ( i.e. neglecting gluon mass corrections) and choosing the 
values^ fi/rriQ = 0.6, 1.0 and 2.0 yields the values of a s {m c ) and a s (m&) reported in 
Table 10. It can be seen that there is poor agreement with the world average values [[|: 

a s (1.5GeV) = 0.357iffi, a s (4.9GeV) = 0.217t°;°or- 



These are calculated using Eqn.(9.5) of reference [0 and correspond to four active quark 
flavours. Matching to the five flavour region where A^ 5 ^ = 2191^3 MeV (corresponding to 
the world average value a s {Mz) = 0.119 ±0.002) is done using Eqn.(9.7) of reference jITJ] 
at a matching scale of 4.3 GeV. For both the J/tp and the T the best agreement is found 
for fi/niQ = 2.0, but the respective deviations are still 4.8a and 3.0a. Fig. 7 shows a 
comparison of the a s (m c ) and a s (rrib) values quoted in Table 10 for ju/mq = 1.0 with the 
present world average value of a s (Q) |§. 



6 Eqns.(7.2) and (7.3) have no real solution for a s (/i) when [i/rriQ = 0.5 
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Figure 9: The effect of gluon mass corrections (m g = 0.721 GeV) on the differential decay 
rate dT/dzdx2 (Eqn.(4-1) for J/tp —> 7gg decays, a) phase space effects only, b) inclusion 
also of longitudinal gluon polarisation states. 
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Figure 10: The effect of gluon mass corrections (m g = 0.721 GeV) on the differential 
decay rate dT/dxidx 2 (Eqn.(7.6) for J/ip — > ggg decays, a) phase space effects only, b) 
inclusion also of longitudinal gluon polarisation states. 
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The gluon mass correction factors are calculated by integrating the differential dis- 
tributions of gluon and photon energies of the decay processes V — > 7gg or V — > ggg 
over the kinematically allowed regions of their respective Dalitz plots. For V — > 7gg it 
is convenient to use the photon spectrum given by integrating Eqn.(4.1) over the gluon 
energies 



1 dY 

r dz 



1 



x_ 



where: 



7T Z 



9 L z 2 
ki(x+/x- 



■[8(27/ - l) 2 (2-4r/ + 7r/ 2 ) 



z 2 (-2 + Ar] + z) 3[ 
+8(2?7 - 1)(5 - 12?7 + 10r/ 2 + 2r?)z 
+2(2/7 - 1)(— 17 + 10t? + 6r] 2 )z 2 + 2(- 



277 + 2r] 2 )z 3 ] 



(1/X--1/X+) 

z 2 (-2 + 477 + z) 2[ 1 ' 



1) 2 (1 + 3t? 2 ) + 4(277 - 1)(3 - 4r/ + 2r/ 2 + 2r7 3 )^ 



+2(7 - 18r7 + lOr/ 2 + 10r] 3 )z 2 + 4(2 + rj)(2rj - l)z 3 + (2 + rj)z 4 



(7.5) 



x-t = 1 — 2r) 



l=Fl/l 



4r/ 



For the decays — > ggg a two dimensional integration is performed over the distribu- 



tion 39 : 



1 dT = 1 1 r I? n _ 3 ?(1 _ 51 15 

r d Xl dx 2 dx 3 " (tt 2 - 9) (xi) 2 (4) 2 (x' 3 ) 2 L 3 /l " 1 8 ' 4 M 

+ (K) 2 + (x' 2 ) 2 + (x 3 ) 2 )(l - U V + 487/ 2 + 25r/ 3 ) 
- 2((x' 1 ) 3 + (^) 3 + (4) 3 )(l-^-3r ? 2 ) 

+ ((x;) 4 + (x' 2 ) 4 + (x' 3 ) 4 )(l + |)] (7.6) 
The allowed phase space region is: 



2 = X!+x 2 + x 3 (7.7) 

2^ < xi < I — 377 (7-8) 

x™" < x 2 < x™ a * (7.9) 

= l-j[l-D( Xl , V )] (7.10) 

C = 1-y [l + D(x u rj)} (7.11) 



£>(^) = ) (7.12) 

V xf 1 — x\ + 77 

The allowed regions of the Dalitz plots for J/t/> — > 7gg decays (open contour) and 
J/t/> — > ggg (cross-hatched contour) for m g = 0.721 GeV are shown in Fig.8a. Similar 
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contours for the corresponding T decays and m g = 1.18 GeV, are shown in Fig. 8b. The 
phase space suppression factors due to gluon mass effects are the ratios of the areas inside 
the contours to the area above the dotted line, corresponding to m g = 0. It is seen that 
the phase space suppression is considerable for T and very large for J/-0 decays. 

The effect of the inclusion of longitudinal polarisation states for the gluons is illustrated 
in Figs. 9 and 10, which show decay rates as a function of photon and gluon energies as 
given by Eqn.(4.1) ( J/^ — > ■jgg) and Eqn.(7.6) (J/ 1 ^ - ► ggg), respectively, for m g = 0.721 
GeV. In Fig. 9a, 10a the longitudinal gluon contributions are suppressed by setting r] = 
except in the equations defining the Dalitz plot boundary, i.e. only the phase space 
limitations due to the non-vanishing value of m g are taken into account. In Fig. 9b and 
10b the complete formulae (4.1) and (7.6), respectively, are used. The most dramatic effect 
of the longitudinal contributions is a strong suppression of the decay rate for J/tp — > 7gg 
in the region z ~ 0.0, X2 — 0.0. As the experimental data analysed here has z > 0.3, 
this has no practical consequences for the present work. Indeed, in the region of small z, 
the dominant mechanism of direct photon production is expected to be the fragmentation 
of light hadrons into photons ||38|| , [f4*8"f, which is not taken into account in the NLO 



QCD calculation of Reference ||14" |. As can be seen in Figs. 9b, 10b the other effect of the 
longitudinal contributions is a modest suppression of the decay rate, near the centre of the 
allowed region of the Dalitz plot, relative to the boundaries. The gluon mass correction 
factors given by integrating Eqn.(7.5) over z or Eqn.(7.6) over x% and X2, are presented in 
Table 11. The rows labelled 'LW use the complete formulae, those labelled 'Phase Space' 
have t] = except in the equations defining the kinematic limits. It can be seen that that 
longitudinal gluon effects are negligible in the correction factors for V — > ggg decays for 
both the J/ip and the T. For V — > 7gg decays these effects increase / 7gg by 30% and 
8%, respectively, for the J/i/j and the T. The errors quoted on the correction factors are 
derived from the total errors on m g given in Sections 5 and 6 above. 

The values of a s (m c ) and a s (rrib) derived from Eqns.(7.2) and (7.3), taking into account 
gluon mass effects according to the values of /ggg//7gg given in Table 11 are presented in 
Table 12 (phase space corrections only) and Table 13 (full gluon mass corrections). In 
each case the values for fi/rriQ = 0.6, 1.0 and 2.0 are presented as well as those given by 
neglecting the 0(a s ) corrections in Eqns.(7.2) and (7.3). For all choices of renormalisation 
scale, the agreement with the world average values is improved as compared to the values 
presented in Table 10, where gluon mass effects are neglected. The best agreement (at 
the 0.3-1.5cr level) is found for /i/mQ = 2.0, though almost equally consistent results 
(deviations of 1.2-1. 7a) are found when the 0(a s ) corrections are neglected. The inclusion 
of longitudinal gluon effects, in the latter case, increases the value of a s (m c ) by 30% 
and a s (mb) by 9%. These shifts are comparable to the uncertainties on a s due to the 
experimental errors on R' v and the gluon mass correction factors. Similar shifts are found 
for yu/wiQ = 2.0 and somewhat smaller ones for /i/mQ = 1.0. The values of a s (m c ) and 
ct s (mb) given by using the full gluon mass correction with /V m Q = 1.0 are compared, 
in Fig. 7, with the world average value of a s (fi), and the values obtained for the same 
renormalisation scale, but without gluon mass corrections. 

Comparison of Tables 12 and 13 with Table 10 and inspection of Fig. 7 shows that that 
the inclusion of effective gluon corrections are essential in order to obtain values of a s (m c ) 
and a s (mb), derived from measurements of R v , that are consistent with the current world 
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average determination of a s (Q). 

Each of the published experiments extracted a value of a s from the measured values 
of R' r . In all cases, good consistency was found with other available measurements of a s . 
This is due to the relatively large errors, at the time, both on the individual measurements 
and on the average value of the other measurements with which they were compared. It 
does not at all contradict the results shown in Table 10, which show instead a poor 
consistency of values derived from the weighted average value of R' r of the three latest 
experiments with the current world average value of a s . The experiments CUSB, ARGUS, 
Crystal Ball and CLEO all used the Brodsky-Lepage-Mackenzie (BLM) scale setting pro- 
cedure [HJ to determine a s at a scale of 0.157M T , i.e. 1.5 GeV. These results, together 



with their weighted mean, are presented in Table 14. The mean value of a s (1.5GeV) of 
0.228 ± 0.019 differs from the current world average value of 0.357 + 0.013 - 0.019 by 4.8 
standard deviations, and is consistent with the results for a s (rrib) given in Table 10. The 
CLE02 experiment used the Principle of Minimal Sensitivity (PMS) to determine: 

a s (M r ) = 0.163 ± 0.002 ± 0.009 ± 0.010 

where the first error is statistical, the second systematic, and the third due to the estimated 
uncertainty of the PMS scale setting procedure. Evolving to the scale mj = 4.9 GeV using 
Eqn.(7.4) gives 

a s (m b ) = 0.190™ 

This differs from the current world average value cited above by only 1.3<r, but also agrees 
within 0.19-l.lcr with the values quoted in Table 10. The latter, however, differ from the 
world average by 3.0-5.7<r. The results of the present analysis, and the combined average 
of those already published in the literature are thus in agreement, and lie 3cr or more 
below the world average. No consistent description is obtained unless the effective gluon 
mass effects are taken into account. 



8 Discussion 



The inclusive photon spectrum for J/ip decays (dashed curve) and T decays (solid 
curve) obtained from the fits performed here to all available experimental data are shown 
in Fig. 11. In both cases, the relativistic corrections and the Photiadis HO QCD correction 
are included. The values 0.721 GeV and 1.18 GeV of m g found in Sections 5 and 6 above, 
are used for the J/ip and T respectively. It is clearly seen that the end-point suppression 
is much more severe for the J/ip than the T. Also shown in Fig. 11 is the RDF spectrum, 
which has been found to describe well all the measurements of the photon spectrum in 
T — > jX except that of CUSB. It is seen to be in good qualitative agreement with the T 
fit curve, but to predict a much harder spectrum than the fit curve for the J/ip. In fact, 
the RDF model, where the gluon mass is perturbatively generated using a low cut-off 
value of 0.45 GeV, predicts that the spectra are very similar in shape for the J/ip and the 
T. This is clearly not the case. 

The very different shapes of the spectra for the J/ip and the T can only be understood 
if the scale introduced into the kinematics of the process by the effective gluon mass is not 
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Figure 11: Inclusive photon spectra including gluon mass effects. Dashed curve: J/if> — > 
7X ; m g = 0.721 GeV ; solid curve: T — > ^X, m g = 1.18 GeV. The dot-dashed curve 
shows the RDF f&$ prediction for T — ■> 7X. 



small in comparison with the mass of the J/^/\ This condition is very well satisfied, since 
the rest mass of the two effective gluons of 2 x 0.721 = 1.44 GeV is 47% of the J/if> mass. 
As previously discussed in MCJHF1 and MCJHF2, the stronger suppression of the end 
point in J /if) decays can be understood as a propagator effect acting on off-shell gluons 
if the genuine gluon mass is ~ 1.0 — 1.5 GeV, i.e. larger than the fitted value of m g for 
the J /if), and similar to that found for the T. Corroborative evidence for this picture is 
provided by the structure of the hadronic final state. The dominance, for massive gluons, 
of the process gg — ► qq over gluon splitting gg — > gggg, ggqq, qqqq leads to a similar 
hadronic final state X in radiative J /if) decays, to the annihilation process e + e~ — > qq at 



the same energy JT3[ consistent with the experimental observations. 



The huge difference observed in the shape of the J /if) and T spectra in Fig. 11 is 
clearly at variance with the principle of Local Parton Hadron Duality (LPHD) where 
parton level pQCD calculations are used (as in the RDF model) down to scales of a few 
hundred MeV ||52|| . Indeed, in Monte Carlo models that give a good detailed description 
of hadronisation effects 0] the cut-off scale of perturbative QCD effects is in the 
range 1-2 GeV, comparable to the effective gluon mass in T decays, and much larger than 
Aqcd- This point will be further discussed below. 

Some remarks are now made on the related work presented in MCJHF1 and MCJHF2. 
In this case fits were performed only to the MARKII data for the J /if), and to those of 
ARGUS and Crystal Ball for the T. In the fit to the J /if), only phase space gluon 
mass corrections were included, with no HO QCD or relativistic corrections. The T fits 
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used phase space gluon mass corrections and the Photiadis HO QCD correction, but no 
relativistic correction. The values obtained for m g of 0.66 ± 0.0l[| GeV and 1.17 ± 0.08 
GeV, respectively, are similar to those 0.721 + 0.016-0.068 GeV, 1.18 + 0.09-0.29 GeV, 
found in the present paper. The larger errors quoted here result from a study of theoretical 
systematics (relativistic corrections, different HO QCD corrections) not done in MCJHF1 
and MCJHF2. Due to a programming error, the resolution functions used in these papers 
had a width that was too large by a factor v^2- This had the effect of destroying the 
sensitivity of the high z part of the J/ip spectrum to the process J/ip — > jrf. A fit with 
an acceptable Confidence Level was then obtained without explictly taking into account 
this decay channel as described in Section 5 above. An important difference between 
the present work and MCJHF2 is an improved understanding of the effect of relativistic 
corrections. As is clear from the discussion at the beginning of Section 2 above, the 
'binding energy' E(p) — tuq introduced by KM must be a positive definite quantity. The 
same conclusion can be reached from simple physical reasoning. In the presence of the 
relativistic correction the heavy quark-antiquark annihilation process occurs over a finite 
spatial region around the origin of the radial wave function, instead of at the origin as 
in the static limit. As the ground state wave function peaks at the origin, relativistic 
corrections must always reduce the decay rate, not increase it. In MCJHF2, following 



the NRQCD [|T^] approach, the corresponding parameter r was taken to be free, to be 
determined from experiment, and was allowed to take positive or negative values. In the 
present paper the relativistic correction parameter < v 2 > is set to the fixed values 0.28 
and 0.09, respectively, for the J/ip and T on the basis of potential model calculations. 
Finally, the NLO QCD correction to the inclusive photon spectrum Jl4| was not available 
when MCJHF2 was written. 

The analysis presented in this paper has neglected possible Colour Octet contribu- 
tions to the radiative decay rates. These have been calculated for T decays at NLO in 
the NRQCD formalism by Maltoni and Petrelli [55|. In the region of interest for the fits 



performed in the present paper, z > 0.3, the corrections to the LO spectrum were found 
to be modest, ~ 10 — 15%. In a more recent study []56 in which octet operators were 



resummed to yield the so-called shape functions |57| a much larger contribution was pre- 



dicted in the near end-point region. However, comparison with the CLE02 data showed 
that the colour octet contribution, with normalisation fixed by the velocity counting rules 
of NRQCD, exceeds the experimental data by between one and two orders of magnitude. 
It may also be remarked that the result of the shape function calculation, in which clus- 
ters of 'non-perturbative' soft gluons are summed, is expected to be drastically affected 
(reduced) if the phase space suppression associated with a gluon mass of ~ lGeV is taken 
into account. In view of the small correction to the shape of the spectrum due even to 
the full colour singlet NLO correction as compared to that resulting from gluon mass 
effects (see Figs. 3-6), the neglect of possible colour octet contributions is not expected 
to modify, in any essential way, the conclusions of this paper. In fact, according to Ref- 
erence |56| the colour octet contributions are expected to strongly enhance the rate in 
the end-point region, whereas what is observed is actually a strong suppression. Indeed, 
both colour octet and colour singlet NLO QCD corrections should be redone, taking into 
account gluon mass effects. Phase space restrictions are expected, in this case, to signifi- 
cantly reduce the NLO corrections to both colour singlet and colour octet contributions, 



7 There is a misprint in MCJHF1, propagated also to MCJHF2, where this error is wrongly quoted as 
0.08 
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Author 


Reference 


Estimation Method 


Gluon Mass 


Parisi, Petronzio 


im 


J/-0 -> 7 X 


800 MeV 


Cornwall 




Various 


500 ± 200 MeV 


Donnachie, Landshoff 




Pomeron parameters 


687-985 MeV 


Hancock, Ross 


m 


Pomeron slope 


800 MeV 


Nikolaev et al. 






Pomeron parameters 


750 MeV 


Spiridonov, Chetyrkin 


o 


n^, (TrGl) 


750 MeV 


Lavelle 





qq — > qq, (TrG 2 ul/ ) 


640 MeV 2 /Q(MeV) 


Kogan, Kovner 





QCD vacuum energy, (TrG 2 ) 


1.46 GeV 


Field 


m 


pQCD at low scales (various) 


1.51^ GeV 


Liu, Wetzel 


m 


lip, (TrGl„) 
Glue ball current, (TrG 2 ^) 


570 MeV 
470 MeV 


Yndurain 


m 




QCD potential 


10~ 10 -20 MeV 


Leinweber et al. 


m 


Lattice Gauge 


1.02 ± 0.10 GeV 


Field 


This paper 


J/-0 -> 7 X 
T -> 7 X 


0.721™ GeV 
1.18±8 gl GeV 



Table 15: Estimates of the value of the gluon mass from the literature. For Donnachie 
and Landshoff, the inverse of the correlation length a is quoted. 



especially for J/ip decays. 

The values of the effective gluon mass determined, in the present paper, from radiative 
J/i/j and T decays are now compared with the results of other studies in the literature of 
gluon mass effects. A number of representative estimates of the gluon mass are presented 
in Table 15. In the following the generic symbol M g will be used for the genuine gluon 
mass, reserving the symbol m g for the 'effective mass' in the sense described in Section 4 
above, determined at tree level in the radiative decays of heavy quarkonia to light hadrons. 



Pioneering work in this field was done by PP 12|. An estimate of ~ 800 MeV for 
the gluon mass was made from the observed softening of the end-point of the inclusive 
photon spectrum in radiative J/ip decays. The work presented in this paper is, in essence, 
a more refined version of the analysis of PP, taking into account experimental resolution 
effects, and the including T decays as well as the best current knowledge on relativistic 
and HO QCD corrections. 

The first extended theoretical discussion of gluon mass effects within QCD was made 
by Cornwall M. Estimates were made of the possible value of a dynamically gener- 
ated gluon mass by several different methods: phenomenological glueball regular isat ion, 
considerations based on the gluon condensate, the glueball spectrum and lattice gauge 
calculations. In fact, many of the gluon mass estimates shown in Table 15 are based on 
the use of the gluon condensate: 

(TrGl) = (0\a s T [G^(x)G^(0)] |0> = N# A (8.1) 



introduced by Shiftman, Vainstein and Zakharov (SVZ) f5"8fl . Here G^ v is the gluon field 
tensor. The gluon mass is related to the inverse of the correlation length, a, of the gluonic 
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vacuum fields. In Reference [58] the phenomenological determination of (TrG 2 ^) by the 



use of QCD Sum Rules is described. 



Donnachie and Landshoff |p9|j , identified the Pomeron trajectory, used to describe 
diffractive scattering, with the QCD two gluon exchange process ||60|1 . They modified 



the perturbative gluon propagator in the long distance region by introducing a finite 
correlation length. The value of the latter was derived from phenomenological Pomeron 
exchange parameters. In Table 15 the reciprocal of the correlation length is equated 
to the gluon mass. Other studies of the sensitivity of the Pomeron parameters to non- 
perturbative modifications of the gluon propagator were made by Hancock and Ross 



and Nickolaev et al. p2| . Again, an effective gluon mass of about 800 MeV was found. 



Spiridonov and Chetyrkin |6^] estimated the gluon mass by calculating power cor- 
rections to the polarisation tensor U e ™(q) of the electromagnetic current of light quarks, 
and identifying them with the gluon condensate term in the Operator Product Expansion 
(OPE) for this quantity derived by SVZ. The calculation was later repeated by Liu and 
Wetzel |p9|| , who obtained a very similar, though not identical, result (see below). 



Lavelle [63], by considering a long-distance modification of the gluon propagator in 
the amplitude for quark-quark scattering, established a relation between a running, dy- 
namically generated, gluon mass and the gluon condensate of SVZ. It was emphasised by 
Lavelle that the derived formula, based on the OPE is only valid in the deep Euclidean 
region. Even so, in a recent paper Mihara and Natale (MN) |65| applied the running 
gluon mass formula of Lavelle to decays of the and T to ggg of 7gg where the gluons 
are either on shell or have time-like virtualities. MN concluded that the average effective 
gluon mass should be smaller for T than for J/ip decays, at variance with the results of 
the fits presented in the present paper. The correction factors due to gluon mass effects: 
/ggg were calculated for the J/ip and T and found to be 0.47 ± 0.30 and 0.94 ± 0.03 
respectively, to be compared with the values found here (see Table 11) of 0.18+0.08-0.02 
and 0.61+0.16-0.06. The smaller gluon mass corrections found by MN are unable to ex- 
plain the large differences in the values of a s (mq) determined from R! v and the the world 
average value of a s (see Table 10 and Fig. 7). 

The gluon mass estimates presented in Table 15 are, with one exception, in the range 



from a few hundred MeV to about 1.5 GeV. The exception is the paper of Yndurain [36] 
which claims that experimental upper bounds in the range from 20 MeV to 10~ 10 MeV 
may be set. These limits are based on considerations of the quantum mechanical potential 
between a quark and an anti-quark. This is assumed to be Coulombic for short distances r 
(r <C (Aqcd)^ 1 ) where Aqcd is the QCD scale parameter, linear ( ~ Kr, K ~ (0.5) GeV 2 ) 
for (Aqcd) 1 < r < m^ 1 and, finally, for r > m~ , to exhibit a Yukawa form: ~ e rrrlg /r. 
The linear portion of the potential gives rise to a barrier of height E cr u ~ Km~ l . Yndurain 
argues that E cr i t may be identified with the highest energy at which unsuccessful searches 
for liberated quarks have been performed. For example, E cr i t ~ 200 GeV leads to the 
limit : 

ro « = 205b = 2 - 5 MeV - 

Other limits are given by applying similar arguments to the absence of proton decay 
into free quarks ( m g < 20 MeV), and the non-observation of free quarks on cosmological 
scales ( m g < 2 x 10~ 10 MeV). The arguments leading to these limits are clearly untenable 
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because of the neglect of the Quantum Field Theory (QFT) aspects of the problem^. In 
fact, when a colour field is stretched between a quark and an anti-quark, the number of 
colour charges is not conserved. The vacuum energy materialises as quark-antiquark pairs, 
which form bound states of light mesons ( 7T, p ,...) This mechanism, as implemented in 
the JETSET Monte Carlo Program |53| is found to give a good description of the observed 
hadron multiplicity in e + e _ — > qq — > hadrons, where the final state is just a colour singlet 
qq pair of the type discussed by Yndurain. This neglect of the QFT aspects of the problem 
renders it unecessary to discuss further the contradiction between Yndurain's upper limits 
and the other gluon mass estimates in Table 15. It is interesting to note, however, that 
the paper of Yndurain is the only one cited on the subject of experimental limits on the 
gluon mass in the current Review of Particle Properties || . 

The estimations of the gluon mass based on the phenomenologically determined value 



of the gluon condensate |58j cited in Table 15 are, with one exception, in the range from 
500-700 MeV. As discussed above, such a value for the gluon mass could not explain 
the much stronger suppression of the end point of the inclusive photon spectrum in J/tp 
decays, as compared to T decays. The exception, that is much more consistent with 
the gluon mass value suggested by the radiative decays, is the calculation of Kogan and 
Kovner ||67|| . This uses an analytical approach in which the vacuum energy of a gauge- 
invariant QCD wave functional is minimised. It leads to the relation, for a pure SU(3) 
Yang-Mills theory^: 

(TrG%) = ^M* (8.2) 

Substituting the value of (TrG 2 ^) of 0.038 GeV 4 from Reference @ gives M g = 1.48 
GeV. This value is quite consistent with the effective gluon masses derived earlier Jl5|, [U| 
and confirmed in the present study. It also agrees well with a previous, independent, 
estimate of the present author, based on a pQCD analysis of several processes with low 
physical scales U68~| , and a recent Lattice Gauge estimate . 



An effective gluon mass has also been introduced in the context of the estimation of 
power-corrections to various hard QCD processes [|7(J . This was already briefly discussed 



in MCJHF2. The leading correction to the mean thrust in e + e~ annihilation into hadrons, 
for example, is found to be ~ M g / ^/s. The gluon mass appears as an intermediate param- 
eter in these calculations, but no explicit values are quoted, and to date, no comparisons 
have been made with other estimates of the effective gluon mass, e.g. lattice calculations. 
Power corrections are an example of 'higher twist' effects in the language of the OPE 
formalism. In fact, it was pointed out in Reference |58j that, as a consequence of the 



dimensionality of the relevant operators, physical systems described by an OPE (as in 
the case of QCD Sum Rules) are expected to have a leading power correction ~ l/Q 4 . 
A OPE description, as used in QCD Sum Rules, is not appropriate to describe higher 
twist effects in the quarkonium radiative decay spectrum discussed in the present paper. 
As can be seen by inspection of Eqns.(7.5) and (7.6) the leading higher twist effects are 



rrig/Q) where Q = M> 



v- 



8 Interestingly enough, Yndurain mentions, near the end of his paper, the screening of the potential 
by quark pair production, but does not draw the conclusion that the existence of such effects invalidate 
his limits on m g derived from a pure quantum mechanical potential. In the real world it is impossible to 
deconfine a quark since the production of low mass qq pairs (mesons) is always energetically favoured. 

9 i.e. effects of quark fields are neglected 



40 



The gluon mass estimate labelled 'pQCD at low scales' in Table 15, is based on phe- 
nomenological arguments (similar to those used later in Reference ||67|), proposed in Ref- 
erence JTIJ and further developed in References |)8|, |7^, |73| . In this approach, effective 



gluon and quark masses are related to the QCD A parameter, that plays the role of an 
infra-red cut-off of the theory. 

It is interesting to note that the approach just mentioned, which relates the A parame- 
ter of pQCD to effective quark and gluon masses, is complementary to the SVZ approach 
of describing the perturbative/non-perturbative interface in terms of QCD Sum Rules. In 
this context it is interesting to consider the relation between ot s , (TrG 2 ^ u ) and M g [ |33T| : 

(TrGD = ^M 4 (3 In |J + f - 12C(3)) (8.3) 

The value M g = 750 MeV quoted in Table 15, was obtained by setting Q = 10 GeV in 
Eqn(8.3), but the equation is expected to be valid for any value of Q in the perturbative 
region. Substituting the values Q = 2.88 GeV, a s (Ql) = a s (0) = 0.27 and M g = 1.0 
GeV, which reproduce well the measured value of a s (M§) p5| , into Eqn(8.3) gives, for 
the gluon condensate, (TrG 2 ^) = 0.094 (GeV) 4 , which may be compared with the average 
value derived by Narison [0 of (TrG 2 ) = 0.071 ± 0.009 (GeV) 4 . Liu and Wetzel m 



derived a formula identical to Eqn(8.3) except that the term 23/2 is replaced by 10. With 
the same parameters as quoted above, the value (TrG 2 ^) = 0.053 (GeV) 4 is obtained. In 
view of the ~ factor of two difference between the original SVZ estimation of the numerical 
value of (TrG^), as compared to that of Narison, the overall consistency of the value of 
(TrG 2 ) derived from the QCD Sum Rule approach, and the pQCD phenomenology of 
Reference [68|, is very satisfactory. 



Thus, although gluon mass effects are only directly observable in processes with pure 
gluonic parton-level final states, such as the radiative heavy quarkonia decays analysed 
in detail in this paper, the existence of a gluon mass of order one GeV is already implicit 
in all pQCD analyses that use, as an infra-red cut-off, the conventional A parameter. It 
is important to stress, however, that as previously pointed out ]F!|, 74], the scale at which 



pQCD is expected to break down is ~ gluon mass ~ 1 GeV, not Aq C d, which is typically 
a factor of five smaller. It is interesting to note that the same conclusion has recently been 
reached in a study of the five-loop QCD /3-function, using Pade-approximant methods pi 



From this point-of-view the success of pQCD, in association with the LPHD hypothesis, in 
describing observed particle multiplicity distributions, as well as their energy dependence, 



using cut-off scales Q as low as 270 or even 150 MeV [p2j] , seems somewhat mysterious 
and may be accidental. As previously mentioned, Monte-Carlo generators that actually 
simulate in detail both the partonic and hadronic phases of the space-time evolution of 
the final state use infra-red cut-off parameters of about 1-2 GeV, of the same order as the 
observed gluon mass. 

Since the main effect of the inclusion of a phenomenological gluon mass on the QCD 
predictions is phase spaee suppression, it is clear that the associated mass must be time- 
like: m 2 > 0. This behaviour is also consistent with an analytical parameterisation of 
the lattice results of Reference |){J . Another recent lattice study using a Coulomb gauge 



gluon propagator [76 has suggested rather a pure imaginary pole mass: 2 (575 ± 124) 



MeV. This work is related to the suggestion of Gribov |77j that colour confinement is 
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due to a long range Coulomb force. In this model physical (transverse) gluons disappear 
form the physical spectrum in the infra-red region. This would appear to be at variance 
with the essential role of massive, time-like, and dominantly transverse, physical gluons 
in the description of inclusive photon spectra that is demonstrated in the present paper. 
Actually, renormalisability remains to be proved for the Coulomb gauge and, as pointed 
out by the authors of Reference [ff6| , except possibly in a confined phase, in quantum 
field theory the square of a particle mass identified with the pole of a propagator must be 
real and positive. The author's opinion is that pure quantum field theory studies of the 
type carried out by Gribov and Zwanziger J7S1 are no better adapted than the quantum 



mechanical potential model of Yndurain pGfl , discussed above, to understand the physical 



mechanism of confinement. This appears to occur via a transition beween partonic and 
hadronic phases of matter, after which colour charges are all confined within hadrons, 
which are mostly bound states of quarks. Clearly quenched lattice calulations of the type 
done in Reference ]76[ are unable to describe such a mechanism. 



It has recently been pointed out that the introduction of a term containing a tachyonic 
gluon mass A with A 2 ~ —0.5 GeV 2 in the OPE of QCD Sum Rules is able to explain some 
long standing puzzles in the related phenomenology J7j|. This approach was justified 
by evidence from lattice gauge calculations for non-perturbative contributions leading 



to a linearly increasing term in the static QCD potential at short distances |3(J. The 
connection of this gluon mass parameter, which gives an economical description, within 
the QCD Sum Rule formalism, of a non perturbative short distance effect with the time- 
like effective mass discussed in the present paper (and all references cited in Table 15), 
is unclear. This latter mass describes rather the long distance behaviour of the gluon 
propagator. It is quite possible that both the time-like and tachyonic gluon masses may 
be appropriate and consistent phenomenological parameters within their different domains 
of applicablity. 



9 Summary and Outlook 



In this paper a phenomenological QCD analysis has been performed using all available 
data on the inclusive photon spectrum in the decays of the J/ip (Mark II Collaboration) 
and the T (CUSB, ARGUS, Crystal Ball and CLE02 Collaborations). The fits performed 
to the shape of the spectra included, for the first time, the combined contributions of 
relativistic corrections, NLO pQCD corrections, and corrections due to the non-vanishing 
gluon mass. For the relativistic correction, fixed values of < v 2 > of 0.28, and 0.09 
were assumed for the J/ip and T repectively, in accordance with recent potential model 
calculations. Both the new, complete, NLO pQCD calculation by Kramer jnj of the T 
photon spectrum and the old resummed calculation of Photiadis [|17| (applicable only in 
the end-point region z ~ 1) were used in the fits. Gluon mass effects were estimated using 



the complete tree-level calculation of V — > 7gg of Liu and Wetzel [39 



For neither the J/ip nor the T was any consistent description of the experimental data 
possible in the absence of gluon mass corrections. In this case, for the Mark II data, no fit 
was obtained with a confidence level greater than 10~ 30 (Table 5). For the T, acceptable 
confidence levels of 0.17 and 0.13 were found for some fits to the Crystal Ball and CLE02 
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data, but not for ARGUS where the best confidence level obtained was 6.7 x 10 3 (Table 
7). However the best confidence level for the combined fit to these three experiments was 
only 4.7 x 10~ 3 . 

Including gluon mass corrections yielded confidence levels of greater than 1% for all 
the fit hypotheses tried (see Tables 6 and 8) and values of the effective gluon mass, m g , 
of: 

m g = 0.7211°;™ l^GeV (J/V) 

(Mark II) 
m g = l.lSSg j£gGeV (T) 

(Mean of ARGUS, Crystal Ball and CLE02) 

The first errors quoted are experimental, and correspond to the 68% error contour of the 
fit, while the second errors are theoretical, reflecting the spread in the best fit values of 
m g resulting from different treatments of relativistic and HO pQCD corrections. It is 
clear from the fit results shown in Tables 6 and 8 that the shape of the photon spectrum 
for both the J /if) and the T is completely dominated by gluon mass effects. Introducing 
the relativistic correction leaves the fitted values of m g almost unchanged, while the shifts 
produced by different HO pQCD corrections are less than, or comparable to, the fit errors. 
The gluon mass corrections are due essentially to phase space restrictions. Including or 
excluding contributions from longitudinal gluon polarisation states changes the fitted 
values of m g by only ~ 5%. 

It can be seen, from the theoretical photon spectra calculated including the effects 
given by the experimentally determined values of m g (Fig. 11), that the suppresssion 
of the end point of the spectrum is much stronger in the case of the J/ip, than that of 
the T. This is in contradiction with the prediction of the QCD parton shower model of 
R.D. Field [R3 that agrees well with the measured T spectrum, but predicts, for the J/ip, 



a much harder spectrum than that actually observed. Such a strong suppression of the 
J/ip end-point, in comparison with that of the T is only possible if the gluon mass is ~ 
1 GeV. All these conclusions are in agreement with those of two earlier, closely related, 



papers [|15|, [16[ , in which it was conjectured, following the work of Parisi and Petronzio [12 
that gluon mass corrections must be much more important than either relativistic of HO 
pQCD corrections in determining the shapes of the inclusive photon spectra. 

The QCD coupling constants a s (1.5GeV) and a s (4.9GeV) were determined from the 
experimental measurements of the branching ratio R' v (Eqn. 7.1) for the J/ip and T 
respectively. Use of R' v has the advantage of being insensitive to relativistic correction 
effects. As shown in Fig. 7, even allowing for a variation of the renormalisation scale in the 
range from 0.6mQ to 2.0mQ, poor agreement (deviations of > 4.8a, 3. Oct for the J/ip, T 
respectively) is found with the current world average value of a s (Q). Using the measured 
m g values to calculate the gluon mass correction factors in the theoretical expressions for 
R' v (Eqns. (7.2) and (7.3)) leads to values of a s that are consistent, albeit within much 
larger errors, (due mainly to the large theoretical error on m g ), with the expected value 
of a s (Q) (see Fig. 7). 

The results obtained in this paper for the effective gluon mass in radiative J/ip and 
T decays are compared, in Table 15, with some other estimates of the gluon mass that 
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have appeared in the literature over the last 20 years. Apart from the upper limits of 
Yndurain ||66|| , that have been critically discussed in the previous Section of this paper, all 
the estimates lie in the range from 500 MeV to 1.5 GeV. The values of m g obtained here 
favour somewhat higher values of the genuine gluon masss M g of > 1 GeV. In particular, 
there is good agreement with the estimates of Kogan and Kovner ||67[| , (minimisation of 
the energy of a pure Yang-Mills QCD wave functional), a previous estimate of the present 
author, using a pQCD analysis of several processes with low physical scales, and the recent 
Lattice Gauge estimate of Leinweber et al. [^] with improved lattice sampling in the far 
infra-red region. It is also mentioned in the previous Section that the conventional value 
for the phenomenological scale parameter A of pQCD of several hundred MeV actually 
implies a gluon mass some five times larger. Since for scales less than M g , A is actually a 
(calculable) scale dependent parameter |ffl| , there is no Landau pole in QCD at the scale 
Q = A, but rather a break-down of pQCD at a much larger infra-red cut-off scale Qo > 1 
GeV. 



A model of confinement in which a purely imaginary gluon mass is introduced [[76 1 
predicts the decoupling of physical gluon states at low physical scales. Such behaviour 
would seem to be inconsistent with the radiative decay data which can only be described 
by the contribution, in the same infra-red region, of physical gluons with a time-like 
effective mass. 

A contribution in the OPE of QCD Sum Rules corresponding to a tachyonic gluon 
mass was recently proposed [f^| to describe some short distance properties of the QCD 



potential suggested by lattice studies [p0[ . This work has no obvious connection with 
the present one where it is shown instead that a gluon with a time-like effective mass 
apparently plays an important role in the long-distance region. 

Since the theoretical uncertainties on the values of the effective gluon mass, m g , deter- 
mined here are much larger than the experimental ones ( the former are 9% (J /if)), 23% 
(T) as compared to 1.3% (J/VO, 5% (T) for the latter) the most urgent need is for im- 
proved theory predictions rather than more precise experimental dataQ. As a first step, 
NLO pQCD calculations should be repeated including gluon mass effects, as has been 
done at tree level by Parisi and Petronzio and Liu and Wetzel. Of even more interest may 
be the calculation of: 

J ftp -> 7qq 

via the exchange of two massive virtual gluons (actually a NNLO process) to test the 
conjectured dominance of this process as already indicated by the structure of the hadronic 
final state. 



As was already evident from previously presented results | I6| , the usefulness of the 
NRQCD approach of Reference |18| will be very limited unless some means is found 
of incorporating the numerically very important gluon mass effects within the formal- 
ism. Another important problem that this approach must address is the possible double 
counting of pQCD and relativistic corrections |73|. 



10 It is interesting to note that the effective gluon mass is determined with a much better relative 
precision for the J/?p than for the T. This is a consequence of the larger size of the gluon mass corrections 
for the 3/tjj resulting in a greater sensitivity to m g . 



44 



Finally it should not be forgotten that the most dramatic effects discussed in this 



paper (in the process J/ip — > 7X) are based on the results of a single experiment [pL3l 
performed more than 20 years ago now. It is clearly important that this remarkable 
experimental result, strongly suggesting that the gluon mass is ~ 1 GeV should be 
confirmed. In a world where half dozen or so b-factories (or their equivalents) exist, or 
are under construction, it is high time that the r charm-quark energy region be revisited 
with modern general-purpose detectors and much higher integrated luminosity than that 
which yielded the data from the Mark II collaboration that has been analysed in this 
paperQ. In view of the less than perfect consistency of the different measurements of 
T — > 7X, an improved measurement of this process, with at least an order of magnitude 
greater statistics, and reduced systematic errors, would be of great interest. This is 
perhaps possible at existing b-factories. 

11 It seems, at the time of writing, that there is now a very good chance that this will soon be done [K2| 
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Appendix 



In the previously published analysis of the and T inclusive photon spectra [16 
resolution effects were simulated by smearing the theoretical distributions with a gaussian 
random number. In order to obtain stable fit results, an unsmeared histogram with a 
very large number of entries was necessary, and the random number throwing had to be 
repeated at each fit iteration, which was both time consuming and resulted in different 
statistical errors in the fitted function at each iteration. Even with ~ 10 6 events in the 
histogram, the determination of the exact position of the \ 2 minimum was quite difficult. 
To avoid these problems, a new, purely analytical smearing technique was used for the 
fits presented in this paper. The method was faster and eliminated all problems related 
to the statistical errors inherent in Monte Carlo methods. 

A theoretical histogram of the inclusive photon spectrum, with bin index J, H0(J), 
is generated with a fine binning, typically NBIN = 1000. The corresponding resolution 
smeared histogram with bin index JS, HS(JS), is the generated according to the following 
algorithm: 

NBIN 

HS(JS) = Yl H0(J)F(J, JS) 
j=i 

where 

(t[z{J)\ 

= 0.0 A(J, JS) > 6, A (J, JS) < -6 



and 



A (J, JS) 



z{J) - z{JS) 



a[z(J)\ 

Here z(J) is the central value of the scaled photon energy in bin J and cr[z(J)] is the 
photon energy resolution at energy E 1 = z(J)My/2. Values of cr[i£ 7 ]/.E y for the different 
experiments analysed here are reported in Table 2. In this way, a gaussian resolution 
smearing of the photon energy by up to ±6<r around the true value is performed. Precise 
results were obtained by normalising the unsmeared histogram to a million events: 

NBIN 

H0(J) = 10 6 

J=l 

but, unlike in the case of a Monte-Carlo simulation, the execution time of the program 
does not depend on this number. The bins of the smeared histogram are grouped to 



correspond to those of the experimental histogram before fitting ]5I |. In all the fits the 
relative normalisation of the experimental and the smeared theoretical histogram was 
allowed to float. 
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